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We introduce a scheme that combines photon-assisted tunneling by a moving optical lattice with strong Hub¬ 
bard interactions, and allows for the quantum simulation of paradigmatic quantum many-body models. We 
show that, in a certain regime, this quantum simulator yields an effective Hubbard Hamiltonian with tunable 
bond-charge interactions, a model studied in the context of strongly-correlated electrons. In a different regime, 
we show how to exploit a correlated destruction of tunneling to explore Nagaoka ferromagnetism at finite Hub¬ 
bard repulsion. By changing the photon-assisted tunneling parameters, we can also obtain a t-J model with 
independently controllable tunneling t, super-exchange interaction /, and even a Heisenberg-Ising anisotropy. 
Hence, the full phase diagram of this paradigmatic model becomes accessible to cold-atom experiments, de¬ 
parting from the region t ^ J allowed by standard single-band Hubbard Hamiltonians in the strong-repulsion 
limit. We finally show that, by generalizing the photon-assisted tunneling scheme, the quantum simulator yields 
models of dynamical Gauge fields, where atoms of a given electronic state dress the tunneling of the atoms with 
a different internal state, leading to Peierls phases that mimic a dynamical magnetic field. 
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I. INTRODUCTION 

Quantum many-body physics studies systems of interacting particles governed by the laws of quantum mechanics. This task 
becomes particularly challenging in a variety of contexts in which the interactions induce strong inter-particle correlations. For 
instance, this strongly-correlated behavior appears in condensed-matter models whenever the system cannot be divided into 
weakly-interacting parts, such that the whole cannot be understood as a sum of its parts and perturbative methods become 
futile [1]. This inherent complexity underlies the abundance of interesting phases of matter that emerge at different scales, but 
also the difficulty in understanding them from an original microscopic model (e.g. high-Tc superconductivity [2]). The same 
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occurs at much higher temperatures and densities, where quarks and gluons interact strongly, and lead to a variety of phases that 
defy our current understanding (e.g. quark matter [3]). In the opposite regime, that of extremely low temperatures and densities, 
ultracold atomic gases trapped by electromagnetic fields are gradually becoming a paradigm of strongly-correlated behavior in 
quantum many-body physics [4]. In contrast to the above condensed-matter and high-energy scenarios, ultracold atoms have 
a unique property: their microscopic properties can be fully characterized and controlled in experiments. This experimental 
control has reached such a status that the dream of exploiting a quantum system to understand the properties of a complex 
quantum many-body model (i.e. a quantum simulator [5]) is already an experimental reality [6]. 

Ultracold gases of neutral atoms can be trapped in periodic optical potentials obtained from the interference of laser beams. 
The dynamics of the atoms in these optical lattices resembles that of tightly-bound electrons in metals, such that this system 
can be considered to be a synthetic solid whose dimensionality and lattice structure can be experimentally tailored, while the 
nuisance of impurities, disorder, or other uncontrolled microscopic degrees of freedom present in real solids, is totally absent. 
Starting from this synthetic solid, it is possible to design a variety of quantum many-body models whose microscopic parameters 
can be experimentally characterized and controlled. For instance, the scattering of atoms leads to a short-range interaction that 
can be tuned all the way from weak to strong repulsion, such that the superfiuid-insulator quantum phase transition of the 
bosonic [7-9] and fermionic [10-12] Hubbard models becomes accessible to experiments. For sufficiently strong repulsion, the 
half-filled Hubbard model leads to a Heisenberg antiferromagnet [13], which yields a playground for quantum magnetism with 
two-component bosonic atoms [14, 15], and the starting point to study high-Tc superconductivity with fermionic ones [16, 17] 
upon controlled doping (i.e. inserting atomic vacancies with respect to the half-filled system). 

In this article, we will combine this strongly-correlated behavior with external periodic drivings to obtain a flexible quantum 
simulator of quantum many-body models. In the context of optical lattices, there is a large body of relevant results regarding 
periodic drivings by modulations of the trapping optical potential. For instance, it is possible to periodically modulate the 
phase of the laser beams forming the optical lattice, as already demonstrated in experiments of chaotic dynamics with cold 
atoms [18, 19]. Another possibility is to modify the detuning of these laser beams linearly in time, usually referred to as 
lattice acceleration, which leads to a linear gradient (i.e. constant force) in the lattice reference frame, and gives rise to Bloch 
oscillations [20]. From this perspective, the previous phase modulation [18, 19] may also be interpreted as a periodic forcing. 
The combination of these two forces permitted probing the Wannier-Stark ladder spectrum [21], and testing the phenomenon of 
coherent destruction of tunneling in the absence of the gradient [22, 23]. Recently [24], a state-dependent coherent destruction of 
tunneling has been demonstrated in optical lattices by using a modulated magnetic field gradient instead of the phase modulation. 
In presence of an energy gradient, one study the phenomenon of photon-assisted tunneling [22, 25]. We shall be particularly 
interested in such photon-as sis ted tunneling effect, whereby the atoms can tunnel in the presence of an energy penalty (i.e. the 
linear gradient) by absorbing photons from the external driving (i.e. the periodic phase modulation). 

Photon-assisted tunneling (PAT) by phase modulation has also turned out to be a useful tool for quantum simulations. The 
dependence of the dressed tunneling on the modulation parameters has been used to drive the system across the superfiuid- 
insulator transition [26-28], and to control the tunneling anisotropy of Bose-Hubbard models in triangular lattices leading to 
magnetic frustration [29, 30]. A subject of research that has received considerable attention recently is the quantum simulation 
of orbital magnetism, whereby the atoms mimic the behavior of electrons in solids subjected to additional magnetic fields [6]. 
Since the atoms are neutral, one must design specific schemes to simulate the effect of artificial/synthetic magnetic fields [32, 33], 
and PAT by phase modulation has also been exploited in this respect (see [31] for a recent review that also covers schemes that 
do not exploit PAT). When the phase modulation leads to an inhomogeneous periodic forcing [34], it is possible to dress the 
tunneling with an effective complex phase. Unfortunately, this simple proposal does not allow for the quantum simulation of 
synthetic magnetic fields [35], and alternative schemes have been considered. For instance, two-tone phase modulations lead to 
synthetic fluxes in arbitrary lattices [36], while the simpler single-tone phase modulations yield staggered fiuxes in certain types 
of lattices [37-39]. 

A different possibility would be to abandon the periodic phase modulation, and investigate other types of drivings that can lead 
to the aforementioned synthetic Gauge fields. Instead of modulating the phase of the optical lattice, one can introduce a periodic 
driving by considering a bi-chromatic deep optical lattice, which can lead to staggered synthetic fiuxes [40] . Alternatively, a 
simple periodic driving by using a pair of slightly detuned and weaker Raman beams (i.e. a shallow moving optical lattice), 
which has been considered in the context of PAT for trapped ions [41], ultracold atoms [42], and generic lattice models [43] 
that can be applied to a variety of contexts. For ultracold atoms, this moving optical lattice yields, in a certain regime, an 
inhomogeneous periodic modulation of the on-site energies of the effective Hubbard model, which can be exploited as a fiexible 
PAT toolbox for quantum simulations of synthetic Gauge fields [41-45]. Here, the atoms tunnel in the presence of an energy 
penalty (i.e. again, a linear gradient) by absorbing photons from the external periodic driving (i.e. this time, the moving optical 
lattice), and acquire a Peierls phase that plays the role of a synthetic Gauge field, and depends upon the wavevectors of the 
Raman beams. 

In this work, we explore a modification of this scheme by considering that the energy penalty can also be caused by the 
on-site Hubbard interactions, yielding a Hubbard blockade that inhibits the tunneling of atoms involving double occupation of 
a lattice site. The combination of this Hubbard blockade with the periodic driving by a moving optical lattice will induce an 
interaction-dependent photon-assisted tunneling. Let us note that the interplay of Hubbard interactions, linear gradients, and 
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phase modulation of the optical lattice, has been shown to be responsible for interaction-shifted resonances in the PAT of Bose- 
Hubbard dimers [46] and chains [47] . Similar effects have been observed experimentally by considering a periodic modulation 
of the intensity of the optical-lattice laser beams [48, 49], rather than the aforementioned phase modulation. This interaction- 
dependent PAT can lead to new schemes to control effective magnetic Hamiltonians [48, 49], or to methods that enhance the 
effects of three-body interactions [50]. We should also mention other proposals that are relevant for the particular subject of our 
work. These concern the engineering of density-dependent tunnelings by either combining laser-assisted schemes with state- 
dependent lattices [51] in the spirit of the original proposal [32], or a periodic modulation of the Hubbard interactions [52-56]. 

In this work, we will show that the interaction-dependent PAT by a moving optical lattice offers a very flexible quantum 
simulator for paradigmatic models of strongly-correlated electrons, and can even allow for the quantum simulation of synthetic 
Gauge fields that are dynamical, in contrast to the static ones mentioned above. As explained below, such synthetic Gauge fields 
evolve under a free Hamiltonian that is not Gauge invariant, and thus depart from the standard Lattice theory approach to Gauge 
theories. 

This article is organized as follows. In Sec. II, we introduce the scheme to implement the interaction-dependent PAT with 
ultracold atoms in optical lattices, and derive a set of effective Hamiltonians that depend on the specific driving, lattice dimen¬ 
sionality, and fermionic/bosonic quantum statistics. The scope of the many-body phenomena that can be studied through these 
effective Hamiltonians is discussed in Sec. III. Finally, we present our conclusions and outlook in Sec. IV. 


II. INTERACTION-DEPENDENT PHOTON-ASSISTED TUNNELING 


In this section, we present a detailed proposal to combine PAT by periodic drivings with strong Hubbard interactions in 
experiments of ultracold alkali atoms in optical lattices. We show that by controlling (i) the atomic interactions by Feshbach 
resonances, and (ii) an additional moving optical lattice, one can exploit an interaction-dependent PAT to delve into interesting 
quantum many-body models that arise in the condensed-matter and high-energy scenarios. 

The starting point is, as customary [4, 57], a trapped atomic gas described in second quantization 
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where we set /i = 1 henceforth. Here, (r), (r) create-annihilate atoms with mass at the position r, and in the electronic 

state I g) corresponding to a particular energy level £a of the atomic groundstate manifold. To remain as general as possible, we 
consider that the components labelled by a may correspond to the states of a bosonic gas, a fermionic one, or a mixture of both, 
which will determine the particular algebraic relations of the creation-annihilation operators. 

We have introduced an optical trapping potential Vot(^) = La^o,a sin^(krcf) + that consists of: (i) A state- 

independent periodic potential, where Vo^a are the ac-Stark shifts of independent pairs of retro-refiected laser beams of wave¬ 
length A =2nlk, which propagate along the axis a G {v,y,z}, and are far detuned with respect to the excited atomic states. To 
obtain state-independent potentials, we assume that the detunings of the laser beams with respect to the excited states are much 
larger than the energy splittings in and that the retro-refiected beams along each axis have parallel linear polarizations [58]. 
To obtain independent potentials along each axis, the pairs of interfering beams must have orthogonal polarizations, or detuned 
frequencies, with respect to other pairs of beams propagating along a different axis. Therefore, it is possible to tune the lattice 
depths Vb,a independently by controlling the beam intensities, which allows to tailor the effective dimensionality of the system. 
(ii) A harmonic trapping caused by a combination of the laser Gaussian profile and the retro-reflection scheme, where 
is the characteristic trapping energy assumed to be sufficiently weak co^ a ^ ^a a — where E'r <7 = k^/lma 

is the so-called recoil energy. 

The final ingredient of the cold-atom Hamiltonian (1) is the ^--wave scattering, which dominates at sufficiently low temper¬ 
atures. This is described by a contact pseudo-potential (r — r') = 47ta(j(jf5{r — characterized by the reduced 

masses = rnamfjf/{rria and the scattering lengths a^yfyi for the collisions of two atoms in the internal state |cr, &). 

Such scattering lengths can be modified experimentally through an external magnetic field via the so-called Feshbach reso¬ 
nances [59]. In the following sections, we show how to exploit an interaction-dependent PAT as new tool to engineer quantum 
many-body Hamiltonians by tuning these scattering lengths appropriately in the presence of a weak moving optical lattice. 


A. Scheme for a periodically-modulated ultracold Fermi gas 

Let us consider a single-species gas of fermionic atoms with two hyperfine states |t) = |F,M), \\) = such that there 

is a unique mass = m^=:m and recoil energy E’r ^ = E’r | =: E’r. We introduce the Wannier basis, (r) =Iiw(r-Ri)/i,<j, 
where w(r — Ri) are the Wannier functions, and /i <7 are the fermionic operators that annihilate an atom of pseudospin a = {t, i} 




4 



® (0\ 



U^ 


n 


tx _ir^ 


\ 

i 


% 


X/2 i 3A'/2 a 


172 A 3A/2 2X 


V 


Vo,. 


Jo 


Jo 


5c/n Jo 


Jo 


1/2 i 3A'/2 2'A 



Figure 1: Scheme of the spin-independent PAT for fermions: (a) Laser scheme corresponding to the static optical lattice formed by the 
retro-reflected beams (red arrows) of frequency 0 ) 2 , and the moving optical lattice formed by a slightly detuned Raman beam of frequency (Oi 
(blue arrow) and the counter-propagating laser beam of frequency of frequency oyi (red arrow). All the laser beams have a linear polarization 
perpendicular to the plane of the figure, as illustrated by the filled circles, (b) Same as before, but considering that the detuned Raman beam 
(blue arrow) forms an angle with respect to the static optical lattice beams (red arrows), (c) Fermionic atoms in two hyperfine states |t) 
(green circles), ||) (orange circles) are trapped at the nodes of a static optical lattice potential (red lines). In the regime of strong 5 -wave 
scattering, atoms tunnel with strength tx between unoccupied sites centered at the energies Conversely, tunneling of one atom to an 

already-occupied site is inhibited by the large energy penalty tx <C U^i (i.e. Hubbard blockade), (d) Moving optical lattice potential (blue 
lines for snapshots of the wave traveling at speed v = Aco/Ak). The tunneling involving doubly-occupied sites can be reactivated when the 
atoms absorb r photons from the moving lattice, providing the required energy to overcome the interaction penalty. (e,f) The PAT contains 
diagonal (e) and off-diagonal (f) correlated events. The diagonal terms correspond to a dressed tunneling within the subspaces of single- or 
doubly-occupied sites, and is controlled by the Bessel function Jq := Jo(t7). The off-diagonal tunneling connects the subspaces of single- 
and doubly-occupied sites, and is controlled by the Bessel function := Jr(t7) and dressed by the tunneling phase Also note that the 

residual dressed interaction is changed with respect to the original bare one 


at the minima of an optical lattice potential Ri labelled by the vector of integers i. We shall consider cubic optical lattices, 
although we note that the scheme detailed below can be directly applied to any other lattice geometry. In this basis, the general 
Hamiltonian (I) can be expressed in terms of the standard Fermi-Hubbard model [II], namely 


= Moc + J^kin + Vint = - II +H.C.) + ^ 

i,c7 i,a <7 ^ i.a 


( 2 ) 


where we have introduced the unit vectors e^, and the notation a = for cr = {t,!}- Here, Ei a = Pa +La 
includes the hyperfine energies and the weak parabolic trapping potential, ta is the tunneling strength of atoms between neigh¬ 
boring potential wells along the a-axis, and Uaa = U^a stands for the on-site interaction due to s-wave scattering, which 
only allows for interactions between fermions of a different state. As customary, we have neglected long-range tunnelings and 
interactions, which requires sufficiently deep optical lattices {Vb,jc: Vqj, Vo £’r. This can be justified considering that 
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while longer range terms are exponentially suppressed with the distance by exp{—^(/?i /A} according to a Gaus¬ 

sian approximation. Let us also note that, for o\^a ^ ^a a’ the harmonic trapping does not modify the tunneling, but simply 
leads to a local term in the Wannier basis that has been incorporated in the local on-site energies £{ (y of the Fermi-Hubbard 
model ( 2 ). 

We consider the limit of very strong repulsion ^ ta, such that the bare tunneling events connecting single-occupied sites 
to doubly-occupied ones are energetically inhibited, as depicted in Fig. 1(c). We shall refer to this tunneling suppression as a 
Hubbard blockade by reminiscence of the Coulomb blockade that inhibits the sequential tunneling of electrons through quantum 
dots. The idea is to overcome this Hubbard blockade via the phenomenon of PAT (i.e. the fermions obtain the required energy 
for tunneling by absorbing photons from an external periodic driving). As shall be shown below, the tunneling of fermions 
between two lattice sites will depend on the density of fermions of the opposite pseudospin populating those sites, which shall 
be exploited to build a quantum simulator. We now discuss two possible periodic drivings that lead to such PAT, and organize 
our presentation by introducing the less demanding schemes first, adding more complexity gradually. 
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1. Two-component fermions in spin-independent moving optical lattices 


(i) One-dimensional scheme: To introduce the main ideas in the simpler setting, let us start by considering a one-dimensional 
(ID) Fermi-Hubbard model obtained from Eq. (2) for {Vbj, Vb,^} ^ Vb,jc, such that only tunneling along the x-axis is relevant 


^FH — ^loc ^kin “b ^int — ^ 

i,a 


E +H-C-) + l'E^aaflaflcfi,afi,a- 


( 4 ) 


As an external periodic driving, we consider a moving optical lattice stemming from a pair of non-copropagating laser beams 
along the v-axis. These beams are slightly detuned with respect to each other (i.e. traveling wave as opposed to the standing wave 
of the static optical lattice. Fig. 1(d)), but again far detuned with respect to the excited states (i.e. Raman beams). Moreover, 
they have the same linear polarization as the laser beams of the static optical lattice to ensure an spin-independent potential [58]. 
Since this moving lattice could induce a spurious tunneling due to the recoil kick imparted by the lasers, we assume that its 
intensity is much weaker Vq ^ (i.e. in the Gaussian approximation). In this regime, the 

effect of the moving lattice is a periodic spin-independent modulation of the trapping frequencies of each potential well 


= £ y COS (AfcX; - Aft)/ + (p) y 


( 5 ) 


where Ak = (ki — k 2 ) • is the wavevector difference, Xf = stands for the minima of the original optical-lattice potential. 
Am = mi — (02 is the detuning of the laser beams, and (p is the relative phase with respect to the static optical lattice. By setting 
Am U^i/r for a positive integer r G Z, the above Hubbard blockade for ^ ta can be overcome through the absorption of 
r photons from the periodic driving (see Fig. 1(d)). To be more precise, as the driving comes from a two-photon ac-Stark shift, 
the process involves absorbing r photons from one laser beam and subsequently emitting them onto the other laser beam. 

To provide explicit expressions for this interaction-dependent PAT, we move to the interaction picture with respect to Uo{t) = 
Sf (exp{i Jq dT(yint +F/niod('Z^))}) 5 such that the fermionic annihilation operators become 


^ Q-'‘Un%o £ J„ (-] f. 

nez 


( 6 ) 


where we have gauged away an irrelevant phase by transforming the fermion operators [60] . The second part of the equality is 
obtained after introducing the important parameter 


r]=yo/Am, (7) 

and using the Jacobi-Anger expansion for first-order Bessel functions J^(z), namely [61]- For simplicity, 

we set XkXi = (| AkA)/ = Ki, which can be achieved with laser beams of the standing and moving lattices of the same wavelength, 
and both propagating along the v-axis. In this configuration, it thus suffices to add a single laser beam detuned with respect to the 
optical-lattice laser beams (see Fig. 1(a)). However, this could be generalized to XkXi = {^XkX)i = Ki/r, which may be relevant 
if the detuned Raman beam does not propagate along the v-axis, but makes some angle with respect to that axis (e.g. r=2 for an 
angle a = K/6, see Fig. 1(b)). 

By substituting the expression (6) in the kinetic Hamiltonian = Uo{t)H^-^ul ((), one finds 

Hkin(r) = -E + H.C.) , (8) 

i,G 

where we have introduced the population difference operator 

= «i+l.CT - (9) 

and a dynamical dressing function 

f(r) = £j„ 

n,m ^ V 2 / 


As announced earlier, the tunneling that connects single-occupied sites to doubly-occupied ones, yielding = ±1, is 

negligible in the absence of the driving Vq = 0 - In this limit, the dressing function is f(() = 1 , such that the dressed tunneling can 
be neglected {tx,a{t)) = ^ 0 in a rotating-wave approximation for tx ^ (see Fig. 1(c)). By switching on the periodic 
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driving Vq 7 ^ 0. this tunneling becomes assisted by the harmonics of the dressing function that are close to resonance with the 
Hubbard interaction, namely for the integers fulfilling = {n — m)Aco (Fig. 1(d)). In particular, by assuming that 

tx, = (U-fi - rA(0) < « rAco, (11) 

we can neglect the majority of tunneling events using a similar rotating-wave argument, except for those that satisfy n = 
+ m. Accordingly, the dressing function becomes simplified 

f(0 = iJm (f) (I) (12) 


We further assume that the laser detuning is chosen in such a way that r is an even integer, and set = 1 for any popu¬ 

lation difference. Making use of the Neumann-Graf addition formula for Bessel functions [61], namely Y.neZ^n{z)Jn^y{z)Q^^^ = 
Jy (2|zsin(0/2) I) we can express the PAT in terms of a single Bessel function 




(13) 


which should be understood in terms of its Taylor series expansion. 

The total time-evolution operator U{t) = can thus be expressed in terms of a time-independent 

Hubbard Hamiltonian of the form (2). However, the dressed tunneling strengths now depend on the density of fermions of the 
opposite pseudospin, and the residual Hubbard interaction depends on the resonance condition in Eq. (11), such that 


Heft = £ ei,aflafi,a - 1 (r? ) flafi+l,a + + ^ £ 


(14) 


As announced in the introduction, the Hubbard-blockaded tunneling becomes activated through a PAT phenomenon, and leads 
to a density-dependent tunneling that can be written as follows 


(^) ~-^0 (t]) (7^1+1,C7 + Jo (t]) ^i,c7^i+l,C7 + Jr (t]) ^i,c7^i+l,(7 + Jr (t]) ^i,c7^i+l,CJ- (15) 

where we have defined the hole number operators = 1 — Ui Q. The first term in Eq. (15) describes the tunneling within the 
subspace of single-occupied sites whereas the second one corresponds to tunneling within the subspace of doubly-occupied 
sites ^ (see Fig. 1(e)). These subspaces can be described as two Hubbard sub-bands centered around Es = 0 and = 5U^^. Fi¬ 
nally, the third and fourth terms stand for tunneling events connecting the single-occupied to the doubly-occupied subspaces (see 
Fig. 1(f)). These four terms can be thus understood as diagonal and off-diagonal tunnelings. We note that a similar classification 
of the tunneling events of the original Hubbard model (4) can be performed by using ^ 

Let us emphasize, however, that the ratio of these diagonal/off-diagonal processes cannot be controlled, which contrasts the PAT 
Hamiltonians (15), where one can adjust the intensity of the moving optical lattice Vq such that the ratio of the Bessel functions 
attains the desired value. This will be crucial to obtain a tunable t-J model with fully controllable parameters in Sec. IIIB. In 
Sec. Ill A, we will use this formulation to connect the effective model to the so-called bond-charge interactions, which leads 
to a quantum simulator of exotic Hubbard models. Moreover, the tunneling of one pseudospin acquires a complex phase that 
depends on the density of the other pseudospin, which will be crucial for the quantum simulation of dynamical Gauge fields in 
Sec. Ill C, when complemented with additional terms that allow us to control each pseudospin independently. 

In order to test the validity of our derivations, we compare numerically the dynamics obtained from the effective Hamilto¬ 
nian (14)-(15), and the periodically driven one (4)-(5) in the simplest setting: a Fermi-Hubbard dimer (see Figs. 2(a,e)). In 
Fig. 2, we explore the real-time dynamics for different configurations of atoms in the initial state, (i) Pauli-blockaded regime: 
Figs. 2(b,c) represent the dynamics of the initial atomic configuration | , ti2) = which does not display Hubbard 

blockade as the tunneling preserves the number of doubly-occupied sites. Nonetheless, the bare tunneling for the spin-up atoms 
(see dashed lines of Fig. 2(b)) is renormalized due to the periodic driving, as shown by the different population dynamics dis¬ 
played by the solid lines (exact) and the symbols (effective). The excellent agreement between the solid lines and the symbols 
proves the validity of our derivations, and the accuracy of the interaction-dependent PAT Hamiltonian (14)-(15). In particular, it 
shows that provided the constrains in Eq. (11) are carefully fulfilled by the system paramaters, terms beyond the rotating-wave 
approximation leading to Eq. (12), within the single-band approximation, do not lead to additional errors departing from the 
desired target Hamiltonian evolution. Regarding the dynamics of the down-spin atoms, we note that these cannot tunnel due to 
the Pauli exclusion principle, as depicted in Fig. 2(c). (ii) Hubbard-blockaded regime: Figs. 2(f,g) represent the dynamics of the 
initial atomic configurations |0i,t|2) = which suffers a Hubbard blockade as the tunneling must change the number 

of doubly-occupied sites. Hence, in the absence of the driving, the atomic tunneling is totally forbidden (see dashed lines of 
Fig. 2(f,g)). By switching on the driving, we observe that the tunneling of both spin-up and spin-down atoms is reactivated, as 
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Figure 2: Interaction-dependent resonant PAT for two-component fermions: Population dynamics of a periodically-driven Fermi-Hubbard 
dimer with parameters tx = 0.\, = 20, and (p = 0 for different density distributions: (a-d) Initial state with one spin-up atom on the right 

well in a background of spin-down atoms | h , t|2) = 1^) ’ initial state with a pair of atoms on the right well |0i, ti2) = 1^) • 

Dashed lines correspond to the tunneling for the un-driven Vq = 0 dimer. Solid lines stand for the resonantly-driven dimer Aft) = t/|^/2 (i.e. 
two-photon assisted tunneling r = 2) with Vq = 3Aft) dictated by the exact Hamiltonian (4)-(5). Symbols stand for the dynamics under the 
effective Hamiltonian (14)-(15) for the same parameters. This criterion of symbols, dashed and solid lines, is kept for other figures, (a, b, c) 
Photon-assisted tunneling for the spin-up atom, and Pauli blockade of the spin-down atoms which cannot tunnel due to the exclusion principle, 
(d) Maximal arrival density of the spin-up atom max{(fzi ^(t)) : 0 < t < 27i/tx}, which displays minima exactly at the zeros of the Bessel 
function Jo(zo,n) = 0 (i.e. coherent destruction of tunneling), (e, f, g) Photon-assisted tunneling for both the spin-up and spin-down atoms, (h) 
Maximal arrival density of the spin-up atom max{(wi |(t)) : 0 < t < 2n/tx}, which displays minima exactly at the zeros of the Bessel function 
h{z2,n) = 0 - In comparison to (d), we observe that the coherent destruction of tunneling depends on the background of spin-down atoms, 
leading to the correlated coherent destruction of tunneling exploited in Sec. III. 


shown by the solid lines (exact) and the symbols (effective), which again show an excellent agreement supporting our analytical 
results. 

Let us now address the phenomenon of correlated destruction of tunneling exploited in Sec. Ill for the quantum simulation 
of strongly-correlated models. According to Eq. (15), the tunneling is dressed by a different Bessel function depending on 
the particle-hole densities, and it can get coherently suppressed when the driving parameter 7] coincides with a zero of the 
corresponding Bessel function. In Fig. 2(d), we observe this effect for the first pair of zeros, rj = zo,n with ^ = 1,2, of the Bessel 
function Jo(zo,n) = 0, which are displayed by the dashed dotted lines. We see how the maximal average population that reaches 
the left site of the Hubbard dimer vanishes when the driving ratio coincides with any of the zeros. In Fig. 2(h), we see that for a 
different particle-hole distribution, the coherent destruction of tunneling takes place at the zeros of a different Bessel function, 
namely rj = Zr,n forn= 1,2 zeros of the Bessel function Jr{zr,n) = 0 for the chosen r = 2. Since the zeros of the two Bessel 
function do not coincide, we can independently suppress the tunneling correlated to a particular particle-hole distribution (i.e. 
correlated destruction of tunneling ), which will be relevant in in Sec. III. 

We have so far presented numerical tests supporting the validity of the resonant PAT, Aft) = U^i/r, such that the residual 
interactions of the dressed Fermi-Hubbard model (14) vanish 5U^i = 0. However, our analytical results show that finite Hubbard 
interactions 5U^i = — rAco) can be achieved by changing the velocity of the moving optical lattice Aft) 7 ^ U^i/ r, which will 

be crucial for several quantum simulations in Sec. III. Let us test this result by numerically integrating an adiabatic evolution 
according to the effective (14)-(15) and periodically-driven (4)-(5) Hamiltonians for a half-filled dimer (Fig. 3(a)). We study the 
evolution of the system (see Fig. 3(b)), for a slow ramp of the tunneling strength tx^tx{l — {Stx)t) with a rate 5tx. Initially, the 
dimer is prepared in the groundstate, which resembles a Fermi sea with the spin-up/down atoms delocalized along the dimer, 
corresponding to the groundstate of the Fermi-Hubbard dimer for 51/^^. After the quench tf^ 1/5^;^;, the dimer should be 
in a spin singlet state corresponding to the groundstate of an antiferromagnetic Heisenberg model that arises for tx 

I'Po) « |FS) = I (4^ + 4) (fl4,+4;) |0) I'Ff) « |HS) = A (f44^-f44) |0), (16) 

In Fig. 3(c), we represent the numerical results for the Heisenberg-singlet fidelity ^s(0 = ^ function of the 

ramp time, and for different ramp rates. We observe that the fidelity approaches ^s(l/ Stx) ~ 1 for the very slow ramps, where 
the adiabatic evolution is expected to be more accurate. Once again, the good agreement between the effective (I4)-(I5) and 
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Figure 3: Adiabatic evolution of an off-resonant PAT dimer at half filling: (a) Scheme of the interaction-dependent PAT for a half-filled 
two-site Fermi-Hubbard model, (b) Scheme for the adiabatic evolution in the half-filled Fermi-Hubbard dimer. At t = 0, the system is initialized 
in the Fermi sea with spin-up/down atoms delocalized along the dimer, which is the groundstate for tx ^ By adiabatically switching off 
the tunneling, the final state for tx <C should be the spin singlet, which corresponds to the groundstate of an antiferromagnetic Heisenberg 
dimer, (c) Singlet fidelity for a periodically-driven Fermi-Hubbard dimer with parameters tx = 0.1, = 20, Vq = 1.84Aft), (p = 0, and 

= O.ltx, with a slow ramp of the tunneling strength tx ^ tx{l — {3tx)t). The symbols correspond to the results given by the effective 
description (14)-(15) for different quench rates 3tx, while the solid lines correspond to the periodically-driven description (4)-(5). 


periodically-driven (4)-(5) Hamiltonians, support our claim that one can study the effects of finite Hubbard interactions, and 
their interplay with the dressed PAT tunneling. 

At this point, it is worth commenting on the effect of higher excited bands that would be present in the optical-lattice setup, 
but are not contained in the single-band approximation implicit to Eq. (2), and the rest of our treatment. The periodic modulation 
may also assist inter-band transition by multi-photon resonances where n/S.(0 = AE, where ^ G Z and AE is the energy gap 
between the lowest and some higher band. To avoid such processes, one must ensure that these resonances are avoided for the 
lowest-lying bands where the number of absorbed photons n can be the lowest. Eventually, this parameter choice may lead to 
a resonance with a much higher band for ^ ^ 1, but provided that AE C the population transferred will be exponentially 
slower than the inter-band tunneling finter-band ^ E(AE/yo)^ ^ E- One must thus make sure that these inter-band population 
transfer is much slower than the time required for the experiment, which becomes essential for the cases that deal with the slower 
super exchange (16). In addition to these possible errors in the simulation, one also has to consider the effect of working in a 
different periodically-modulated picture, which can lead to micro-motion contributions at the driving frequency that can alter 
the experimental measurements [45]. 

Before moving to the PAT in higher dimensions, let us mention that we could gain additional fiexibility in the scheme by 
introducing an additional linear gradient, which may come from a lattice acceleration, an external electric field, or a magnetic- 
field gradient. If we tune the gradient such that it coincides with the on-site interaction, we can generalize Eq. (14) by substituting 

iv) ^ ) ~ {'u) (^ ) ^0 (^ ) "E ^2r (17) 

According to this expression, the off-diagonal tunnelings connecting doubly- to single-occupied sites (i.e. fourth term in Eq. (17) 
represented in the upper panel of Eig. 4(a)), and single- to doubly-occupied sites (i.e. third term in Eq. (17) represented in the 
upper panel of Eig. 4(c)), depend on different Bessel functions. This leads to a two-tone beating in the tunneling dynamics, as 
shown in Eigs. 4(b), (d), which also serve as tests of the validity of our analytical derivations. 

(ii) Higher-dimensional scheme: The scheme presented above can be directly generalized beyond ID. The static optical 
lattice should be modified such that it allows for tunneling along two (Vq,^ ^ {Vb,x 5 Vbj} ^ Er) or three ({Vb,x 5 Vbj, Vb,z} ^ ^r) 
directions. As can be observed from Eqs. (6)-(10), to assist the tunneling along a given direction, it is crucial that the periodic 
modulation (5) has a phase that varies along that particular direction. Therefore, we would need to include additional moving 
optical lattices that propagate along the remaining axes, dressing the corresponding tunneling along two a = {x,y}, or three 
a = {x,y,z} directions. One may consider adding one independent detuned laser beam per axis, paralleling the construction 
of the one-dimensional case. Otherwise, one could simply tilt the laser beam of the one-dimensional case, such that it has a 
non-vanishing projection propagating along each axis. The former scheme would lead to independent moving lattices along 
each axis whose intensity and frequency can be tuned separately, whereas the latter would lead to a non-separable moving lattice 
that dresses all the different tunnellings with the same intensity and frequency, albeit one could play with the propagation angle. 

Eor simplicity, we consider the first situation, such that the periodic driving is 

i,G oc ^ 


( 18 ) 
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Figure 4: Interaction-dependent resonant PAT for two-component fermions in a gradient: Population dynamics of a periodically-driven 
Fermi-Hubbard dimer with parameters tx = 0.1, = 20, Vq = 3.5Aft), (p = 0, and subjected to an additional gradient A = (a) Scheme 

for PAT for states with double occupancies |0i,ti2) = diagram), and | tii 5 02) = /i^^/l^lO) (lower diagram), showing two 

different tunneling dressings that lead to the two-tone beating displayed in (b) for the initial state I'Pq) = |0i,ti2)- As usual, the solid lines 
correspond to the numerical solution of the exact Hamiltonian (4)-(5) in the presence of an additional gradient, while the symbols stand for 
the numerical solution of the effective Hamiltonian (14) with the modified tunnelings (17) due to the presence of the gradient, (c) Scheme for 
PAT for states with single occupancies | ii,t2) = (upper diagram), and | ti 5 i 2 ) = /i|/l^|9) (lower diagram), showing two different 

tunneling dressings that yield the two-tone beating displayed in (d) for the initial state I'Pq) = \ iiAi)- 


where we have introduced the labeling indexes i = {ixjy) for 2D, and i = {ixJyJz) f^^* ^D. As before, we have assumed that 
for 2D (Vo^jc C Vb,x, and ^ and for 3D (Vq^x ^ Vb,jc, <C Vbj, and yo,z ^ Vb,^), such that the moving lattices do not 
modify the bare tunneling and only lead to a periodic modulation of the on-site energies. Each of these moving lattices assists 
the tunneling along a given direction, and does not interfere with the tunnelings along the remaining axes. Accordingly, the 
interaction-dependent PAT is a direct generalization of (14), which requires a parameter regime 


5/7^4, = (C/ti - faAcOa) <C U^i « r:,A(0^ = ryAcOy = r^Aco^, (19) 

and yields the following effective Hamiltonian 

- II (^a) + H-C.) + (20) 


1,<7 


where A^i+ea ,<7 = '^i+ea ,<7 “ the bare tunnelings are approximated by Eq. (3), and the dressed tunnelings depend on 

^raAm+ecc,^ (^la) ^i,a^i+ea,a + Jo (t]a) ^i,a^i+ea,a + fra (^a) ^i,(T^i+ea,C7 + fra (t]a) ^i,(T^i+ea,<75 (21) 

where we have introduced rja = Vq ^/AcOa. It is interesting to note that controlling the intensity difference of each moving 
optical lattice, we can tune the spatial anisotropy of the dressed tunnellings. The possibility of generalizing to 2D is especially 
interesting in the context of the t-J model, and its connection to high-Tc cuprate superconductors, as outlined in Sec. IIIB. 

We have thus seen that the interaction-dependent PAT with a moving optical lattice leads to effective Hubbard models of any 
dimensionality with dressed tunnelings that are density dependent. In the following section, we will show that by considering 
a state-dependent moving optical lattice, the PAT scheme becomes more flexible, which will allow us to target other quantum 
many-body models, in particular dynamical Gauge flelds. 


2. Two-component fermions in spin-dependent moving optical lattices 

(i) One-dimensional scheme: Let us once again start with the less demanding case of ID. We note that the far-detuned moving 
optical lattice can become state-dependent if the laser-beam polarizations are not collinear (see Eig. 5 (a,b)). This occurs even for 
detunings that are larger than the Zeeman and hyperflne splittings, as far as they do not exceed the flne-structure splitting [58]. 
Eor fermionic alkali atoms, this could turn to be incompatible with reaching ultracold temperatures, as the flne-structure splitting 
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Figure 5: Scheme of the spin-dependent PAT for fermions: (a) Laser scheme similar to Fig. 1, but with the linear polarization of the 
detuned Raman beam (blue arrow) rotated by an angle Op with respect to the static lattice laser beams. This gives rise to a spin-dependent 
moving optical lattice, (b) Same as before, but considering that the detuned Raman beam (blue arrow) forms an angle with respect to the static 
optical lattice beams (red arrows), (c, d) The moving optical-lattice potential for each pseudospin (green lines for |t), and orange lines for 
ID) reactivates the tunneling involving doubly-occupied sites. The PAT contains diagonal terms (e) corresponding to a spin-dependent dressed 
tunneling within the subspaces of single- or doubly-occupied sites controlled by the Bessel function Jq^cj := hi'Hcj)- The off-diagonal terms (f) 
contain a spin-dependent tunneling connecting the subspaces of single- and doubly-occupied sites, which are controlled by the Bessel function 
Jr,cj •= Jr(t7cj) and dressed by the spin-dependent tunneling phase 


is rather small, and the residual photon scattering may become appreciable [33]. However, we stress that the moving lattice is by 
construction much weaker than the static spin-independent one. In fact, we can reduce the residual photon scattering by orders 
of magnitude by lowering the intensity of the moving-lattice laser beams, as far as their detuning is simultaneously lowered, 
such that the ratio rj controlling the PAT (15) remains constant. We will thus assume that a conservative state-dependent moving 
optical lattice can be realized without increasing the photon scattering and heating the ultracold atomic gas. 

In this case, we can generalize the driving (5) by including a state-dependent periodic modulation of the on-site energies 

= Y, “D COS {AkXi - Am + (22) 

i,a ^ 

where the spin-dependent amplitude must again fulfill Vo^a ^ and (pa stands for a phase difference with respect to the static 
lattice that is generally state dependent [62, 63] (see Figs. 5 (c,d)). Going back to the spin-independent scheme (5), the new 
modulation (22) can be achieved by rotating the polarization of the laser beam that is slightly detuned with respect to the static- 
lattice lasers [62]. In this case, the spin-dependent driving amplitudes can be tuned by controlling such an angle, or instead 
the direction of propagation of the laser beam with respect to the quantization axis [63]. Another possibility would be to resolve 
the hyperfine structure, such that one can exploit selection rules in the ac-Stark shifts. In fact, for pseudospins corresponding to 
the maximally-polarized Zeeman sublevels, it is possible to obtain optical lattices that selectively address a single pseudospin 
(i.e. Vb,t = / 0)^ or vice versa, as realized in ion-trap experiments [64]. This leads to a spin-dependent driving where the 

wavevector Aka, detuning AcOa, intensity Vo^a^ and relative phase (pa can all be controlled independently for each pseudospin 

^mod(0 ~ COS [Aka^i AcOat (Pa^^i a‘ (^3) 

i,G ^ 

Paralleling the previous section, we will consider equal wavelengths of the static and moving optical lattices, such that Aka^i = 
(^Ak<jA)/ = although we remark again that the scheme also works for other propagation angles. Once the new periodic 
drivings (22)-(23) have been discussed, we can address the interaction-dependent PAT they give rise to. We shall use Eq. (23), as 
the results also encompass those related to the driving (22). By reproducing the steps that lead to the effective Hamiltonian (14) 
for the spin-independent driving, we find a parameter regime analogous to Eq. (11), namely 


5U^l = {U^i - VaAcOa) < ^ 

and a new dressing function of the tunneling that becomes spin-dependent, namely 

fa{t) = (^) ( y) 


(24) 


(25) 
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Figure 6: Spin-dependent coherent destruction of tunneling for two-component fermions: Population dynamics of a periodically-driven 
Fermi-Hubbard dimer with parameters tjc = 0.1, = 20, and subjected to a spin-dependent moving lattice with (p^ = (p^ = 0, = Aft)^ = 

U^l/2. (a) Scheme for the spin-dependent PAT for the state |0i,t|2)^ showing that the spin-down atoms can be coherently frozen, as shown 
in (b) for Vot = 3Aft), and = 5.13Aft) such that 12(5.13) = 0. As usual, the solid lines correspond to the numerical solution of the 
exact Hamiltonian (4) with the spin-dependent periodic modulation (23), while the symbols stand for the numerical solution of the effective 
Hamiltonian (26) with the modified tunnelings (27). (c) Scheme for the spin-dependent PAT for the state |0i, showing that the spin-down 
atoms can be coherently frozen, as shown in (d) for Vot ~ 5.13Aft) such that 12(5.13) = 0, and Vot ~ 3Aft). 


where rfa = /AcDa. In this case, the detunings are chosen such that % is an even integer for both pseudospins, such that we 
can thus set = i Using the Neuman-Graff addition formula once again, we find that 


i,G i,G ^ i,G 


(26) 


Remarkably, we find that the amplitude of the density-dependent tunneling can be controlled independently for each pseu¬ 
dospin 


^raArif^l ^ {^g) — Jo (^(j) + Jo {^g) ^i,G^i^l,G + (^a) (^c) 5 (27) 


and that the tunneling phase of one pseudospin depends on the density of the other pseudospin, which will be crucial for the 
quantum simulation of dynamical Gauge fields in Sec. Ill C. In order to benchmark these predictions, we study numerically 
a spin-dependent coherent destruction of tunneling in a Fermi-Hubbard dimer subjected to the spin-dependent moving optical 
lattice. According to Eq. (27), the dressed tunneling of a doubly-occupied half-filled dimer (see Fig. 6(a)) depends on the spin 
of the atom, such that the spin-up atoms tunneling is controlled by the Bessel function J^^(t]^), whereas the spin-down atoms 
tunneling depends on the Bessel function J^^(t 7 |). Therefore, by controlling the intensities of the the spin-dependent moving 
lattice, the dressed tunneling of the spin-down atoms can be coherently destructed Jr;(^|) = 0 , while the spin-up atoms hop 
freely in the lattice Jr|(^^) 7 ^ 0 (see Figs. 6(a)-(b)). Conversely, we can coherently freeze the spin-up atoms = 0, while 

the spin-down atoms hop freely in the lattice (r]^) 7 ^ 0 (see Figs. 6(c)-(d)). Let us emphasize the excellent agreement between 
our analytical description (symbols), and the exact dynamics of the periodic Hamiltonian (solid lines). 

(ii) Two-dimensional scheme: The scheme presented above can be generalized beyond ID. For the sake of concreteness, and 
for its particular interest in connection to the t-J model in Sec. Ill B, and the dynamical Gauge fields in Sec. Ill C, we will restrict 
to 2D (Vb,z ^ {Vb,x, Vbj} ^ F’r). The idea is to consider spin-dependent moving optical lattices along the x and y axes 


^mod(0 ~ 2 COS {Ak^~ ^^G,ai + ^G,a) 

i,G a ^ 


(28) 


such that the relative phases of the moving lattices fulfill (pa,x > 0, but (pa,y = 0. In this case, and after following the same steps 
as above in an analogous parameter regime 

ixdyi ~ ~ ^ G,a^^G,a) ^ ~ ^ G,a^^G,a — ^ < 7 ', a'^^< 7 ', a': CC (29) 

one can derive the following effective Hamiltonian 


//eff = - I 

i,(7 


i+ex,<7 


■ iy^rajArii^ 


SU^ini,tnu, ( 30 ) 
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Figure 7: Density-dependent Aharonov-Bohm interference for two-component fermions: Population dynamics of a periodically-driven 
Fermi-Hubbard tetramer with parameters tx = ty = Q.l, = 20, and subjected to a spin-dependent moving lattice with (py^a = 0, and 
(Px,a = ^/2, and ^(Oa^o = Va,(7. (a) Scheme for the PAT for the state | til 502503504 )^ where the frozen spin-down atom induces a 

TT-flux in the tunneling of the spin-up atom, leading to the Aharonov-Bohm destructive interference in (b). (c,d) Populations for Vot = 

Voi = 5.13Aft) such that 12(5.13) = 0, and the spin-down atom is frozen. As usual, the solid lines correspond to the numerical solution of the 
exact Hamiltonian (4) with the spin-dependent periodic modulation (23), while the symbols stand for the numerical solution of the effective 
Hamiltonian (30). Due to the interference, the spin-up atom cannot reach the opposite corner = 0. (e-h) Same as above, but for the 

state I til, 125 Is 514)5 where the spin-down atoms cannot tunnel due to the Pauli exclusion principle. For this density background, the tunneling 
phase of the spin-up atoms vanishes since there is no spin-down density gradient. Hence, this initial state does not lead to Aharonov-Bohm 
destructive interference, and the spin-up atom does indeed reach the opposite corner. 


where when atoms tunnel along the x-axis, they acquire a dynamical phase that depends 

on the density of the other pseudospin, whereas they experience a vanishing phase when tunneling along the y-axis. This will be 
equivalent to the so-called Landau Gauge in Sec. Ill C, which is accompanied by a non-vanishing dynamical Wilson loop. 

So far, all of our numerical tests have been independent of the phase of the moving optical lattices, as we have only addressed 
a Fermi-Hubbard dimer. By considering the simplest 2D case, a Fermi-Hubbard tetramer forming a square plaquette, we can al¬ 
ready test numerically the predicted effect of the moving lattice phase, which according to Eq. (30), induces a density-dependent 
Peierls phase in the tunneling. To observe the effects of such a Peierls phase, we shall first exploit the above spin-dependent 
destruction of tunneling, or the Pauli exclusion principle, to freeze the dynamics of the spin-down atoms (see Figs. 6(a-b)). Then, 
the immobile spin-down atoms yield a density background that modifies the tunneling phase of the spin-up atoms. We explore 
this possibility in Fig. 7 for two different density distributions of the spin-down atoms, which lead to the presence/absence of 
an Aharonov-Bohm destructive interference in the tunneling dynamics of the spin-up atom. Besides confirming the validity of 
our analytical description (30) in a transparent scenario, let us note that by lifting the coherent destruction of tunneling, and 
allowing the spin-down atoms to tunnel, the Peierls phase will acquire its own dynamics, which will be crucial for the quantum 
simulation of dynamical Gauge fields in Sec. Ill C. At this point, the reader may skip the following sections, and move directly 
to the quantum simulator applications dealing with fermionic models in Sec. III. 


B. Scheme for a periodically-modulated ultracold Bose gas 


Let us now turn our attention to a single-species gas of bosonic atoms, and consider again two hyperfine states |t) = 
|F,M), ID = such that we have a unique mass =: m and recoil energy = £’r| =: F’r. In the Wannier 

basis, 'Pa(r) = Li w(r — rP)^i,a 5 the general Hamiltonian (1) can be expressed in terms of the bosonic operators as a two- 
component Bose-Hubbard model [8], namely 


Hbh - Hioc +^kin + Lint - + 


1,<7 


i,C7 Ot 


i a,a' 


(31) 


where the Hamiltonian parameters coincide with the fermionic ones (3) under similar approximations. However, due to the 
different statistics, s-wave scattering among pairs of atoms with the same electronic state are now allowed (i.e. D 0), 
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which gives more freedom for the interaction-dependent PAT. We shall discuss different regimes of interest for the quantum 
simulation of bosonic quantum many-body models, which can be achieved by tuning the Feshbach resonances. 


1. Two-component hardcore bosons in spin-independent moving optical lattices 


Let us start from the ID Bose-Hubbard model obtained from Eq. (49) for {Vbj^Vb Vb,jc, such that only tunneling along 
the v-axis is relevant. In the hard-core limit, double occupancies of bosons with the same pseudospin are energetically forbidden 
(i.e. ^ tx). In this limit, we can project out all states with sites occupied by more than one boson of the 

same pseudospin provided that the filling is <1- By using the corresponding projector we can map the bosonic 

creation-annihilation operators onto an 5u(2) spin algebra 










— I 1 — • ^/(T- 


(32) 


In such a hardcore limit, the Bose-Hubbard model (49) only contains the on-site Hubbard interactions for two bosons of opposite 
pseudospin, namely 


^hBH — ^loc ^kin 3“ Vint • — ■ 








H.c)+Wu,^bibyb.ib,^^. 


(33) 


Analogously, we must also project the spin-independent periodic modulation due to the moving optical lattice which, in the same 
regime as discussed for fermions (5), yields 

^mod(0 := ^ ^ COS (AkXi - Acot + (p)n-^ (34) 

i,G ^ 


We can now proceed in analogy to the fermionic gas, as the difference between the fermionic operators and the hardcore-boson 
ones does not change any of the steps of the derivation. We move to the interaction picture with respect to the projected 
driving and Hubbard interactions f/o(0 = ^ /o^'z:(Vint + ^mod('Z ^))})5 and follow the same steps as in the fermionic case 
to express the time-evolution operator IJ{t) = f/J (^)e“^^^* in terms of 


Hcf( = '£ei,abl^h,a 

i,G 


■ (^x^rAHi^ 
i,G 






i,(T’ 


(35) 


l,G 


where the dressed tunnelling strengths and phases are again density dependent. Given the one-to-one correspondence with the 
fermionic scheme, the numerical results to test the validity of our derivation by comparing the full and effective dynamics would 
not add anything different from Fig. 2, and we shall not include them here. 

Since the hardcore-boson and the fermionic number operators have the same algebraic properties, we can rewrite this density- 
dependent tunnelling strength in complete analogy to the fermionic case 

^ (t]) = Jo (t]) + Jo (t]) + Jr (t]) + Jr (t]) (36) 

where the hardcore hole operator is hi^a = 1 — hi^a. Therefore, the dressed tunnelings for hardcore bosons can be described 
pictorially by Figs. l(e,f), which distinguish events that preserve/modify the double occupancy of the tunneling sites. 

Let us also note that the fermionic schemes for spin-dependent drivings in Sec. II A, and the generalization to higher dimen¬ 
sions, equally applies to the bosonic gas in this hardcore limit. Although we shall focus on the fermionic applications of the 
quantum simulator in Sec. Ill, we emphasize that all the quantum many-body models discussed there have a hardcore-boson 
counterpart, including the dynamical Gauge fields. 


2. Two-component softcore bosons in spin-independent moving optical lattices 

The objective of this section is to relax the hardcore constraint ^ ^ tx, which forbids double occupancies of 

bosons with the same pseudospin. Let us, however, start by understanding the PAT of a single-component bosonic gas described 
by Eq. (49), but restricted to a single pseudospin (e.g. a =t). For notational convenience, we drop the pseudospin index, such 
that the Hamiltonian corresponds to the standard Bose-Hubbard model 

^BH — ^loc ^kin 3“ Vint ~ H-H.C.^ + !)• 


(37) 
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Figure 8: Scheme of the spin-independent PAT for bosons: (a) Bosonic atoms in a single hyperfine state |t) (green circles) are trapped at 
the nodes of a static optical lattice potential (red lines). Tunneling of one atom to an already-occupied site is inhibited by the Hubbard blockade 
tx <^U. (b) Moving optical lattice potential (blue lines) reactivates the tunneling by providing the required energy for the atoms to overcome 
the interaction penalty. (c,d) The PAT can be divided into events that conserve the Hubbard interaction energy (c) controlled by the Bessel 
function Jp := Jo(t7), and those that change it (d), and are controlled by various Bessel functions := J^(t7 ) and the corresponding tunneling 
phases Also note that the residual dressed interaction dU is changed with respect to the bare original one U. 


According to our discussion of the Hubbard blockade tx <C U, the tunnelling of bosons that changes the parity of the occupation 
number is energetically forbidden (see Fig. 8(a)). As customary, we activate this tunnelling by means of the periodic modulation 

^mod(0 ~ , (38) 

i 

given by a moving optical lattice acting on the bosonic atoms (see Fig. 8(b)). Due to the different Hubbard interaction, which 
now involves a single pseudospin, one cannot use the expression in Eq. (6). Instead, we consider the interaction picture of a 
bond operator which, again up to an irrelevant Gauge transformation, can be shown to be 

where Uo{t) = ^ (exp{i/odT(yint + ^mod(T))}) is the interaction-picture unitary. After defining the bosonic population dif¬ 
ference operator Ant^i = — Ui, and using the Jacobi-Auger expansion for each of the time-dependent exponentials, the 

expression of the kinetic energy, = f/o(0^kin^^o (O’ becomes 

//kin(f) = -£ (tx{t)b]b.^^+H.c.'], t,{t) = (40) 

i,a 

with exactly the same modulation function as in Eq. (10). 

We can then proceed by following the same steps as for the fermionic case, to find a parameter regime 

tx, 5U = {U- rAco) rAco, (41) 


and finally obtain the effective interaction-dependent PAT Hamiltonian 


J/eff = Y^eibjb, -E +H.C. 


V ( 1 "i 


(42) 


Here, we observe that the dressed tunnelling depends on the density of bosons of the same pseudospin at the sites involved in the 
tunnelling event (see Eigs. 8(c,d)). Since the boson number per lattice site is not restricted anymore by the hardcore constraint, 
we cannot express it as the simple polynomial (36) quadratic in the density operators, but rather as a highly non-linear term. 
Using the orthogonal projectors onto subspaces with a well-defined difference number of bosons this non-linearity 

becomes apparent 


3r(An/_^l-Kl) (t?) ^ (t]) 


A^b+l 


E ^rf{ri)Y[ 

e'=o ejte 


/2-/'2 


(43) 
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Figure 9: Interaction-dependent PAT for single-component softcore bosons: Population dynamics of a periodically-driven Bose-Hubbard 
dimer with tx = 0.l,U = 20, and (p = 0 for different density distributions. (a,b) Initial state with one atom on the right well |0,1). The dashed 
lines correspond to the tunneling for the undriven Vq = 0 dimer, while the resonant driving Aft) = t//2 (i.e. two-photon assisted tunneling 
r = 2) with Vb = 1.5Aft) corresponds to the numerical solution of the exact Hamiltonian (37)-(38) (solid lines), and the effective one (42)-(43) 
(symbols). We use the same criterion dashed-line/solid-line/symbols in all the figures. (c,d) Same as above for an initial state with two atoms 
on the right well and one atom on the left well |1,2). (e,f) Same as above for an initial state |0,2). (g,h) Off-resonant effects on the PAT for 
the initial state |0,2) by modifying the detuning 5U/tx E {0,0.5,1,1.5,2,2.5} (dark-to-bright coloring implies increasing the detuning 5U). 


where A^b is the total number of bosons loaded in the optical lattice, and we have again assumed that r is an even integer. In 
analogy to the studies for the phase-modulation driving [46, 47], we observe that there will be interaction-shifted resonances 
that correspond to the zeros of the Bessel functions for different density backgrounds. In fact, using our formalism, one could 
derive similar analytic expressions for a phase-modulation driving [46, 47]. For the intensity-modulated lattices of the recent 
experiments [48, 49], the situation is simpler as there can only be one particular occupation that is resonant with the periodic 
modulation of the tunneling matrix element, and no Bessel functions arise. The possibility of engineering Bose-Hubbard models 
with a density-dependent tunneling strength and phase by our PAT scheme can be considered as an alternative to the proposals 
based on periodic modulations of the v-wave scattering length [53, 54]. Let us finally note that it is possible to generalize this 
scheme to higher dimensions, paralleling the fermionic case Sec. II A. 

In order to test the validity of our derivations, we have numerically compared the time-evolution predicted by either the 
effective Hamiltonian (42)-(43), or the periodically driven one (37)-(38) for a Bose-Hubbard dimer in the regime of interaction- 
dependent PAT. In Fig. 9, we explore the real-time dynamics for different configurations of atoms in the initial state: (i) Non- 
blockaded regime: Figs. 9(a-b) represent the dynamics of the initial atomic configurations |0,1) = Zp^lvac), which does not 
display Hubbard blockade as it consists of a single atom. Yet, the bare tunneling (see dashed lines of Fig. 9(b)) is renormalized 
due to the periodic driving, as shown by the different population dynamics displayed by the solid lines (exact) and the symbols 
(effective). The same renormalization occurs for the configuration |1,2) in Figs. 9(c-d), which is not blockaded as the overall 
occupation parity is conserved in the tunneling process. On top of the dressing of the tunneling, we observe a bosonic enhance¬ 
ment which leads to a doubled tunneling rate with respect to Fig. 9(b). (ii) Hubbard-blockaded regime: Figs. 9(e-h) represent the 
dynamics of the initial atomic configuration |0,2), which suffers a Hubbard blockade as the tunneling must change the overall 
occupation parity. Hence, in the absence of the driving, the atomic tunneling is totally forbidden (see dashed lines of Fig. 9(f)). 
By switching on the driving, we observe that the tunneling is reactivated, as shown by the solid lines (exact) and the symbols 
(effective). So far, all these simulations correspond to the resonant PAT, where the parameter regime (41) is achieved for 5f/ = 0. 
In Figs. 9(g-h), we explore the off-resonant case 5U > 0, and the possibility of describing the driving detuning as a residual 
Hubbard interaction. The agreement between the solid lines (exact) and the symbols (effective) in Fig. 9(h), shows that this is 
indeed the case. We observe that, as 5U is increased, the periodic exchange of particles is gradually inhibited, as one would 
expect since the single- and doubly-occupied Hubbard bands become more and more separated in energy. 

As a further numerical proof of the consistency of our effective description, let us explore the phenomenon of coherent 
destruction of tunneling. According to Eq. (43), the effective tunneling is dressed by a different Bessel function depending on 
the densities of the bosonic atoms. For instance, the PAT tunneling of Fig. 10(a) is controlled by Jo(t]), whereas the tunneling of 
Fig. 10(c) is controlled by 1 ^( 17 ). Therefore, whenever the driving parameter 7] coincides with a zero of the corresponding Bessel 
function, the tunneling should get coherently suppressed. In Fig. 10(b), we observe this effect at T] = zo^n forn = 1,2,3 zeros of 
the Bessel function Jo(zo,n) = 0, which are displayed by the dashed dotted lines. We see how the maximal average population 
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Figure 10: Density-dependent coherent destruction of tunneling for single-component softcore bosons: Maximum population reaching 
the left well msLx{{ni {t) :0 <t < Titx} for a periodically-driven Bose-Hubbard dimer with parameters tx = 0.l,U = 20, for a resonant drying 
Aft) = U jl (i.e. two-photon assisted tunneling r = 2) with (jp = 0, as a function of the driving ratio rj for different density distributions, (a-b) 
Initial state |0,1), which displays minima in msLx{{ni{t) :0 <t < n/tx} exactly at the zeros of the Bessel function Jo(zo,n) = 0- (a-b) Initial 
state |0,2), which displays minima in max{(wi {t) :0 <t < n/tx} exactly at the zeros of the Bessel function J2(z2,n) = 0- 


that reaches the left site of the Hubbard dimer vanishes when the driving ratio coincides with any of the zeros. In Fig. 10(d), we 
see that for a different atomic density distribution, the coherent destruction of tunneling takes place at the zeros of a different 
Bessel function, namely 7] = Zr,n for ^ = 1,2,3 fulfilling Jr(zr,^) = 0 for the chosen r = 2. 

Once the interaction-dependent PAT for the single-pseudospin bosons has been understood, and its validity has been checked 
numerically, we can turn our attention to the situation of two-pseudospin bosons without the hardcore constraint 


^BH — ^loc ^kin Pint — 




H.c. 


i a,a' 


fb].b: 

1,0 l,G l,G 1,0 


(44) 


We shall be interested in the regime of tx, where double occupancy of bosons of the same (different) pseudospin 

is allowed (forbidden). In this regime, as the intra-spin interactions do not blockade the tunneling, we only include the inter-spin 
interactions in the interaction picture Uo{t) = ^ (exp{i Jq dT(yint,tf+^mod('r))}), where a spin-independent moving optical 
lattice is applied to both pseudospins 


^^mod(0 = E ^ ~ 

i,G ^ 

One can then see that, in the following parameter regime 


(45) 


tx, I/;;, 5U^I = - rAco) < ^ rAft), 


(46) 


the dressed tunnelings will only depend on the density of atoms of the opposite pseudospin, as occurs for the hardcore bosons 
or the fermions. Therefore, the effective Hamiltonian becomes 


Heff = -E ,, +H.C.) 

i,G i,G 


^EE 

i a,o' 




(47) 


where we have introduced the residual interactions and For (p = 0, we get an 

exotic Bose-Hubbard model with the analogue of the fermionic bond-charge interactions, whereby the dressed tunneling of 
one pseudospin depends on all the possible density backgrounds of the other pseudospin through the corresponding Bessel 
functions. The difference with respect to the fermionic bond-charge interactions (15) is the highly non-linear function of the 
bosonic densities, namely 


Nb 

i'=0 


Nb 


e=o 


f - 
£^-£'2 


(48) 


It is also worth commenting that, had we set = Un ^ tx in the blockade conditions, the dressed tunneling would de¬ 
pend on the density of both pseudospins ^ ^ Jka«,+i ^+An,+ 15 + 1 ) (j?)e which might 
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Figure 11: Interaction-dependent PAT for two-component softcore bosons: Population dynamics of a periodically-driven two-component 
Bose-Hubbard dimer with parameters tx = 0.1, = 20, = 0, and (p = 0, for different density distributions. (a,b,c) Initial state 

with one spin-up atom on the right well, and one spin-down atom on the left well 11^, 1|). The dashed lines correspond to the tunneling for the 
undriven Vq = 0 dimer, which shows the Hubbard blockade as a consequence of tx <C 11^^. Switching on the off-resonant driving Aft) 

(i.e. two-photon assisted tunneling r = 2) with detunings 5U^^/tx G {0,0.5,0.75,1} (dark-to-bright transition implies increasing the detuning 
SU^^), we observe that the tunneling for the spin-up atoms (b) (and spin-down atoms (c)) is reactivated, and depends on the residual Hubbard 
repulsion dU^^. Here, the solid lines correspond to the numerical solution of the exact dynamics (44)-(45), while the symbols represent the 
effective dynamics (47)-(48), for a periodic driving of intensity Vq = 3Aft). (d,e,f) Same as above, but for an initial state |1|,2|) and resonant 
driving = 0. Due to the bosonic enhancement of the dressed tunneling of the spin-up atoms (d), the tunneling of each spin-species takes 
place at a different rate. As a consequence, the dynamical Hubbard blockade leads to a more complex tunneling dynamics. 


also be interesting regarding exotic Bose-Hubbard models. Let us note that, once again, this effective description matches 
perfectly the dynamics of the driven two-component Bose-Hubbard dimer (see Fig. 11). 

Once again, we could generalize to higher dimensions, or to spin-dependent drivings, and study the strongly-correlated models 
that arise. However, the properties particular for the restricted number of particles of fermions and hardcore bosons, which allow 
for instance to build a quantum simulator of dynamical Gauge fields (see Sec. Ill), cannot be generalized to the soft-core boson 
case. Let us finally note that, by considering the hardcore-boson limit on Eq. (47), one recovers the previous result (35). 


C. Scheme for a periodically-modulated ultracold Fermi-Bose mixture 


Let us now turn our attention to a mixture of bosonic and fermionic atoms, and consider that the two hyperfine states of 
the original formulation (1) correspond to the fermionic and bosonic atoms |t) = |E,M) =: |f), |i) = |F',M') =: |b), respec¬ 
tively. In the Wannier basis, 'L^(r) = w(r — rP)/i for the fermions, and 'L^(r) = Ei w(r — rP)^?! for the bosons, the general 
Hamiltonian (1) can be expressed in terms of the fermionic f\ and bosonic hi operators [65] as a Bose-Fermi Hubbard model 


//bFH = -£ +H-C.) + £ 

i i,a i 

+ - L (fa.f/iVi+e„ + H.c.) + £t/b44Vi h ■ 


(49) 


Here, the Hamiltonian parameters now depend on the fermionic/bosonic nature of the atoms, such as the tunnelings 


ta,b = b I —^ 
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which differ due to the different masses m{,rrn, through the corresponding recoil energies £R f,£R b. The Hubbard interactions 


Ubb — \j ^^tbb^Rjb I p3 
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( 51 ) 
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not only differ by the mass ratio, but also by the different scattering lengths for a boson-fermion (2bf, and a boson-boson (2bb 
collision. Once again, the fermion-fermion collisions are forbidden by the Pauli exclusion principle. We shall discuss different 
regimes of interest for the quantum simulation of Bose-Fermi quantum many-body models, which can be achieved by tuning the 
Feshbach resonances, and controlling these scattering lengths appropriately. 

Let us focus on the ID case. The moving optical lattices can be made dependent on the bosonic/fermionic atomic species, as 
the corresponding atoms usually have a very different atomic level structure. Therefore, in analogy to Eq. (23), we consider the 
following periodic driving 

^^mod (f ) = L ^ - AcOft + <Pf) f ^ ^ cos {AkbXi - A(Obt + %) n,. , (52) 


where the amplitudes Vo,f/b wave vectors A^f/b, frequencies AcOf/b, and relative phases (Pf/b, depend on the particular atomic 
species, and can be controlled experimentally. These periodic drivings will assist the tunneling against an energy penalty given 
by the Bose-Fermi interaction, which is very large 


— (^bf ^bAtOb) t/bf ~ ^bAtOb — TfAcOf, 


(53) 


where we have introduced the residual Bose-Fermi interactions 5f/bf. and two integers rb,rf G Z that determine how many 
photons are absorbed from the periodic driving to overcome the energy penalty, and assist the tunneling. 

In analogy with the two-component softcore bosons, only the interspecies Hubbard interactions inhibits the tunneling, such 
that the interaction picture must he Uo{t) = ^ Jo d'z:(yint,bf + ^mod('Z ^))})5 leading to the time-evolution operator U{t) = 
Uq (r)e“^^^* with the following effective Bose-Fermi Hubbard Hamiltonian 


= I e^bb]b, - £ (77b) 7.;^ b,^, + h.c.) + W Ubbb]b]b, b, 

i i i 

+H.C.) +\Y,^Ubfb]fJ f, b,. 


(54) 


We thus observe that the bosonic/fermionic dressed tunnelings, both the amplitude and phase, depend on the fermionic/bosonic 
densities. Moreover, they can be independently tuned by controlling the parameters of the bosonic/fermionic moving lattice, 
which will be interesting for the quantum simulator of dynamical Gauge fields (see Sec. III). We also note that the generalization 
to 2D yields a Bose-Fermi analogue of the effective Hamiltonian (30). 


III. APPLICATIONS OF THE PHOTON-ASSISTED TUNNELING QUANTUM SIMULATOR 

We have introduced a scheme to implement an interaction-dependent PAT with trapped ultracold atoms. In this way, we 
have obtained a toolbox with different effective Hamiltonians that depend on the quantum statistics, dimensionality, and spin- 
dependent/spin-independent nature of the periodic driving used to assist the tunnelling. In this section, we will discuss how such 
a toolbox can be exploited for the quantum simulation of interesting problems in condensed matter and high-energy physics. 
Instead of focusing on a particular application, we have decided to stress the wide scope of the proposed quantum simulator by 
describing a variety of interesting quantum many-body models. All these problems share a common feature, they are described 
by complex quantum many-body models, either on the lattice or in the continuum, which still present open questions that defy 
the capabilities of existing numerical methods. We will describe the context of the particular models that can be targeted with the 
quantum simulator, and try to discuss the essence of the phenomena that they try to capture. Moreover, we will highlight what 
we believe are open questions of the models that have been studied in detail over the years, and point to future work required to 
understand the models that have not been explored in such a detail. 


A. Bond-charge interactions and correlated electron-hole tunnelings 

In the original derivation of the Hubbard model for electrons in transition metals [10], it was shown that a variety of longer- 
range terms also arise in the Hamiltonian. In addition to the tunnelling and the on-site density-density coupling. Coulomb 
interactions also yield nearest-neighbor terms that can be described as the repulsive interaction between: (i) charges localized at 
neighboring ions V, (ii) charges localized at an ion and a neighboring bond/linkX, and (in) charges localized at two neighboring 
bonds/links W. These bond-charge interactions can be responsible for a host of interesting effects in the context of charge density 
waves [66], ferromagnetism in itinerant electron systems [67], or alternative mechanisms of superconductivity based on electron 
holes [68], and other interesting effects [69]. However, it has been argued that the required values of A, W with respect to f/,y 
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to observe such effects are not likely to be realized in standard solid-state materials [70]. On the other hand, the possibility of 
controlling these terms in the synthetic solids offered by ultracold atoms in optical lattices has recently raised some interest in 
the community [71]. 

From the perspective of ultracold trapped atoms, one can evaluate these bond-charge terms for the contact interaction in 
Eq. (1). In the Wannier basis, one can introduce the bond operators, flof\ a’ account for the state-dependent 

density located at the bond (i,j). Then, the bond-charge terms modify the standard Hubbard Hamiltonian (2) by introducing 


AN^ = 


EE 

i,a O’ 




^'^i,(7'^i+ea,a+ ^i+ea,a) + 
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,i+eoc i,i+eoc 


(55) 


where we can express the different interaction strengths in the Gaussian approximation as follows 


= \j 






= \l -ka^iER 


vb^Vo,,yo,,V^" 




Wa=Va. (56) 


From these expressions, we observe that the ratios ValU^WalU, and XajU are vanishingly small for deep optical lattices [72]. 
An interesting possibility to reach regimes {Va^Wa^Xa} ^ U, where the bond-charge interactions can lead to new phases of 
matter, is to consider ultracold dipolar gases [71] or, as we show in this section, to exploit the PAT toolbox. 

At this point, we introduce an additional modification of the standard Hubbard model (2), whereby the tunnelling is not only 
modified by the density at separate sites (i.e. bond-charge interaction Xa in Eq. (55)), but also by the density-density correlations 


i,a (y 


(57) 


Here, Xa is the strength of this tunnelling, which cannot be obtained from any two-body interaction in the Wannier basis. In fact 
this term would require rather exotic three-body interactions, which cannot be realized even with ultracold dipolar gases [71]. 
These terms are interesting in a condensed-matter context, where they appear after reducing models with hybridized bands to 
single-band Hubbard models, as occurs in intermediate-valence solids [74], and the high-Tc cuprates [75]. As shown below, the 
PAT toolbox can control all these terms in the effective cold-atom Hamiltonian. 


1. Correlated electron-hole tunnelings: bond-ordered waves and triplet pairing 


In this section, we focus on a quantum simulator of the Fermi-Hubbard model with tunable ratios of tajU^XajU^XalU, 
which may lead to very interesting many-body effects. Such a Fermi-Hubbard model F/fh = ^loc + ^int can be rewritten 

in terms of asymmetric tunnelings that are correlated to the electron/hole occupation 


^corr 

^kin 


EE + ?eh(«i.o'^i+e„,CT + hi;ani+ta,a) + 4“«i.CT«ie„,CT) flcfi,+ta,c + ^-C., 

i,a o 


(58) 


where we have introduced tunnelings in a hole-hole background = —ta, in an electron-hole background t^^ = —ta ^Xa, and 
in an electron-electron background t^^ = —ta +^a + '^Xa. 

Let us focus for simplicity on the ID limit of Eq. (58). The effective cold fermion Hamiltonian (14)-(15), which is ob¬ 
tained through the PAT by a spin-independent moving lattice, already contains these correlated particle-hole tunnelings. For the 
parameter regime fulfilling (11) with r = 2, and setting the moving lattice relative phase to (p = 0, one finds 


4 = -f^Jo(r?), 4 =-r;,J2(J?), t^^ = -txJo(v)- (59) 

Accordingly, the correlated tunneling asymmetry can be tuned by modifying the intensity of the moving optical lattice Vq, or its 
frequency Am, such that 7] = Vo/Aco is varied. This effect could also be achieved through a periodic modulation of the ^--wave 
scattering length as considered in other recent schemes [55, 56]. 

The ID Fermi-Hubbard model with the asymmetric correlated tunnelling strengths (59) hosts a variety of quantum many-body 
phases, contrasting the situation of the standard ID Fermi-Hubbard model, where only insulating and Luttinger-liquid phases 
occur. We now comment on the phases that could be explored with cold atoms given the constraints imposed by the specific 
tunnelling rates (59). This model was initially studied for t^^ = t^^ < t^^, where a groundstate with spin-density-wave order can 
be found for sufficiently strong repulsion at half-filling [80]. Later on, it was realized that for t^^ = t^^ > t^^, a new density wave 
where charge alternates on the bonds (i.e. bond-ordered wave) can be stabilized for not too large interactions [81]. Interestingly, 
it has been recently shown that, by modifying the filling factor, the phase diagram of the model is much richer even for vanishing 
Hubbard interactions [56]. For instance, a superconducting phase with unconventional triplet pairing was identified. In addition 
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to the interest of exploring these phases with cold atoms using the above PAT scheme (14)-(15), this could serve to benchmark 
the proposed quantum simulator, since these results rest on very accurate and efficient analytical and numerical methods that 
exist in ID. Moreover, the quantum simulator could be used to study the fate of the predicted phases for different chemical 
potentials [56], and the appearance of new ones, as the Hubbard repulsion is switched on. 

Once the quantum simulator has been verified, it would be very interesting to consider the 2D Fermi-Hubbard model with 
correlated tunnelings. Although the phase diagram is to the best of our knowledge mostly unknown, the results on the ID model 
suggests that it can host a variety of new phases with respect to the standard 2D Fermi-Hubbard model, the understanding of 
which defies analytical and numerical methods. We believe that the possibility of finding interesting phases of matter, even 
above the stringent temperatures to observe magnetic ordering in the 2D Fermi-Hubbard model, is certainly worth exploring. 

To introduce the topic of the following subsection, we emphasize that it is not possible to set Xa =0 without making Xa =0 
simultaneously, as would be required to study solely the effects of bond-charge interactions (55). This also occurs for the 
schemes based on a periodic modulation of the 5 '-wave scattering length [55, 56]. In the following subsection, we shall show that 
this becomes possible by introducing an additional linear gradient in our scheme. 


2. Bond-charge interactions: hole superconductivity and rj -pairing 


In this section, we focus on a quantum simulator of the Fermi-Hubbard model with bond-charge interactions, whose impor¬ 
tance can be controlled by tuning the ratios of ta/U^Xa/U. For simplicity, we focus on the ID case, and consider the PAT by 
a spin-independent moving lattice (14) in the presence of an additional linear gradient, which leads to the dressed tunneling in 
Eq. (17). For the sake of concreteness, we consider a parameter regime fulfilling Eq. (11) for r = 2, and set (p = 0. By adjusting 
the intensity of the moving lattice to Vq = 1.56A(0 (i.e. = 1.56), such that terms like (57) in the effective Hamiltonian (14) 

vanish X = 0 , we obtain 


ffeff = U^i,a + SUaai)fUi,a 

i,<7 






I,a 


(60) 


where t = tj2 and the bond-charge interaction X = r;^:(Jo (t]^) — J 2 (^★)) alternates between neighboring sites. Other possible 
values of the moving-lattice intensity fulfilling = 0 are T]^ = {4.89,8.29,11.53,...}, and correspond to solutions of the 
equation Jq (t]^) + J 4 (t]^) — 2 J 2 (t]^) = 0. Moreover, one can also change the moving-lattice detuning such that r = {4, 6 , 8 ,...}, 
and look for the solutions of Jq (t]^) + J 2 r (t]^) — 2 J^ (r]^) =0. This will yield different values of r]^, and allow for the tunability 
of the ratio of X/t, and the signs of X^t. This quantum simulator can explore the phenomenon of hole superconductivity [ 68 ], 
where the bond-charge interaction 0 < X ^ r can be responsible for a superconducting phase even in the presence of a repulsive 
interaction 5Uyi > 0, as has been predicted using a mean-field approximation for 2D [76]. Since our effective Hamiltonian (60) 
can be generalized to higher dimensions following our results in Sec. II, the quantum simulator could test the correctness of such 
mean-field predictions [ 68 ]. From a broader perspective, the quantum simulator can explore the phase diagram of the model in 
different regimes, such as X > t, and study the effects of the bond-charge alternation in Eq. (60) 

Since the residual Hubbard interactions in (60) depend on the detuning of the photon-assisted scheme, one can also study the 
attractive case < 0. It has been shown [77] that a bond-charge interaction X = t stabilizes an T] -pairing groundstate [78] for 
finite attraction, which forX = 0 only occurs for strictly infinite interactions 5Uyi —> 00 [ 79 ]. Although our quantum simulator 
cannot reach the exact condition ofX = t (e.g. for = 1.56, Xjt ^ 0.94), it has been argued that relaxing these conditions may 
still host the T] -pairing groundstate, even if the exact methods of [77] cannot be applied any longer. It would be very interesting 
to explore this possibility with the quantum simulator, addressing the role of the finite residual gradient SUao in Eq. (60), the 
bond-charge alternation, and the possibility of achieving a regime X > t that is not feasible for the standard Hubbard model. 

Before closing this section, we note that this can also be addressed with hardcore bosons, following our results in Sec. II. 


B. High-Tc superconductivity and itinerant ferromagnetism 


To introduce the topic of this section, let us consider again the Eermi-Hubbard model (2) with additional bond-charge 
terms (55). The interaction between charges localized in neighboring bonds Wa leads to a couple of terms: (/) a pair tun¬ 
nelling for fermions of opposite pseudospin AH^i = 2 Wa/i^|/j^^/i+ea,t/i+ea,| + H.c., and (ii) a direct exchange interaction 

~ +H.C. between neighboring pseudospin excitations, where we have introduced the spin- 1/2 operators 






^(4/i.t-A/u)' 


(61) 


Regarding the direct exchange, the possibility of observing magnetic ordering of localized pseudospins in the cold-atom sce¬ 
nario would require a direct exchange = 4Wa that is much larger than the attainable values (56). The same occurs for 
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transition metals, where the direct exchange cannot explain the appearance of magnetic ordering. As first pointed out by R W. 
Anderson [13], magnetic ordering can also arise as a consequence of the strong Hubbard interactions that prevent conduction. 
Moreover, in the presence of doping, the interplay of this tendency towards magnetic ordering with the correlated dynamics of 
holes is one of the proposed mechanisms that could explain high-Tc superconductivity [16]. 

In this section, we show that the effective Hamiltonian (20) for the PAT scheme with (jO = 0 may open a new route to study 
models of itinerant quantum magnetism and high-Tc superconductivity, provided that the limit of large repulsive interactions 
is considered. In this regime, the subspaces of single-occupied ^ and doubly-occupied ^ lattice sites become well- 
separated Hubbard sub-bands. The kinetic energy in Eq. (20) contains terms that act within each of these sub-bands 

Ko = {ria)flafi+ea,(y) ^+e „,0 + EE (-tah (ria) fl^fi+ea,<y) «i+e„,a +H-C., (62) 

i,a C7 ^ ^ i,a a ^ 

and terms that connect the two sub-bands, such as the operator ^ ^ ^ expressed as 

i,a <7 i,a (7 ^ 

or the operator ^ ^ which can be expressed as 

^d^s = (~^aha (^a)^^^i+ea,a (~^a^ra iVa) 

i,a C7 ^ i,a a ^ 

We shall use this formulation to propose a quantum simulator that explores Nagaoka ferromagnetism even in the limit of finite 
repulsive interactions, and the full phase diagram of the t-J model in ID and, more interestingly, in 2D. 


7. Itinerant ferromagnetism: correlated destruction of tunnelling and Nagaoka ferromagnetism 

For a strong Coulomb repulsion, electrons in undoped transition metals lower their energy by displaying antiferromagnetic 
ordering. In the standard Fermi-Hubbard model, the origin of this effect can be traced back to the so-called super-exchange 
interaction, whereby an antiferromagnetic spin pattern allows for virtual electron tunnelling between neighboring sites that 
lowers the kinetic energy [13]. The situation is utterly different for ferromagnetism, where the reliability of initial mean-field 
predictions of large ferromagnetic regions in the phase diagram is highly questionable [82]. One of the few rigorous results on 
the existence of ferromagnetism in the Fermi-Hubbard model is due to Y. Nagaoka [83], who showed that a single hole in a large 
class of half-filled Hubbard models can lead to a fully-polarized ferromagnetic groundstate for infinite repulsion. 

The stability of this Nagaoka ferromagnet for different regimes has been a topic of recurrent interest in the literature. For 
two holes [84], the ferromagnet is no longer the groundstate. Nonetheless, for finite hole densities, the fully-polarized Nagaoka 
ferromagnet can be stable up to a critical hole doping [85]. Despite some initial discrepancy regarding its full polarization [ 86 ], 
more recent numerical results based on quantum Monte Carlo [87] and density-matrix renormalization group [ 88 ] agree with 
the above scenario [85]. Since the task of doping a transition metal with exactly one hole seems quite daunting, these results 
are crucial for an experimental realization of the Nagaoka ferromagnet. Another obstacle for the realization of a Nagaoka 
ferromagnet is the requirement of infinite repulsion. In a cold-atom context, alternatives exploiting a long-range double-exchange 
interaction in a two-band Hubbard model [89], or a large spin-imbalance in an optical lattice with a ladder structure [90] have 
been considered. We show below that our PAT scheme allows to access the physics of the infinite-repulsion Hubbard model, and 
thus Nagaoka ferromagnetism, even for finite Hubbard interactions. 

Fet us consider the effective kinetic energy obtained for the 2D or 3D scheme (20). By looking at the tunnelings in Eqs. (62)- 
(64), one notices that by modifying the intensity of the moving optical lattice such that = 0, while Jo(t 7 ^) 7 ^ 0, we can 

exactly cancel the terms that do not preserve the parity in the occupation number, namely Eqs. (63)-(64). For concreteness, 
we consider a parameter regime fulfilling ( 11 ) for r = 2 , such that the moving lattices have the same intensity and detuning, 
and T]^ = Vq ^/AcOa = 5.135. Accordingly, only the parity-conserving tunnelling (62) is preserved. This effect is similar to the 
so-called coherent destruction of tunnelling [91], where the tunnelling of electrons subjected to a periodically-modulated force 
is totally suppressed for certain parameters of the force. In our case, the suppressed tunnelling is fully correlated to a particular 
electron-hole background (see Fig. 2), such that the effect can be understood as a correlated destruction of tunnelling. 

For large but finite Hubbard repulsion the Hubbard sub-bands are separated in energies and no term in the effective 

Hubbard Hamiltonian can connect them. By controlling the atomic filling factor such that + < 1, all the dynamics 

takes place within the single-occupation subspace, and is controlled by 

^eff = =^8^^ (“^«/il(7/i+ea,(7 

i,a <7 
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where ta = and we have introduced the Gutzwiller projector onto the single-occupied sub-band = 

Interestingly enough, the Hamiltonian (65) corresponds to the infinitely-repulsive Hubbard model supporting Nagaoka ferro¬ 
magnetism. However, in our scheme, only a finite repulsion is required to allow for the adiabatic loading of the lower Hubbard 
sub-band. By varying the filling factor, the stability of the ferromagnetic phase and the full phase diagram can be explored 
experimentally. In particular, it can serve as a benchmark of the phase diagram presented [ 88 ], which is based on extrapolating 
numerical results for ladders with increasing number of legs, and has predicted an intermediate phase-separation region between 
the fully-polarized Nagaoka ferromagnet and the paramagnetic phase. 


2. High-Tc superconductivity: tunable t-J and t-XXZ models 


In the previous section, we have considered an alternative route to access the physics of the limit of infinite repulsion in the 
hole-doped Fermi-Hubbard model (i.e. (^i,^ + < 1). However, a large but finite repulsion can also lead to very interesting 

physics. In this regime, the competition of the projected kinetic energy (65) with the antiferromagnetic super-exchange [13] 
leads to the so-called t-J model 

Hu = Hu = EE (-^'^Vi+ea.cT + H.C.) +E/(5i , (66) 

i,a <7 ^ ^ i,a 


where t is the tunnelling within the single-occupied subspace, and / > 0 is the strength of the antiferromagnetic super-exchange 
interaction. In the case of square lattices, this Hamiltonian (66) has been considered as the canonical effective model in the 
theory of the high-Tc cuprates by part of the scientific community [16, 92, 93]. In this context, the t-J model arises after mapping 
a more microscopic three-band Hubbard model of the cuprates [94] onto a single-band one [75, 95]. By using the microscopic 
parameters of the three-band model [96], one finds that J/t ^ 0.3 is the typical regime that can be realized in these materials. For 
the Fermi-Hubbard model with ultracold atoms (2), one can obtain an effective t-J model in the limit of very strong repulsion 
t Then, one finds /= A'f which cannot attain values larger than /< 0.47 since the tunnelling must be at least 

t < O.lUyi to allow for the perturbative process underlying the super-exchange. Therefore, the standard Fermi-Hubbard model is 
almost at the verge of the regime of importance for the hole-doped cuprates J/t ^0.3 [97]. At this point, it should be mentioned 
that reaching the low temperatures required to observe the effect of the antiferromagnetic super exchange at equilibrium is a great 
challenge that has only been achieved recently [98]. However, from a broader perspective, the rest of the rich phase diagram of 
the t-J model cannot be explored with these experiments. A possibility to attain tunability over these parameters, while at the 
same time controlling and homogeneous atomic doping from half-filling, would be to consider composite-fermion quasiparticles 
whose tunneling corresponds to a correlated tunneling of a boson-fermion pair in a Bose-Fermi mixtures [99]. We show below 
that our PAT scheme leads to a standard fermionic t-J model where the ratio J/t can attain any desired value by controlling the 
intensity of the moving optical lattice. In this way, the full phase diagram of the t-J model may become accessible to cold-atom 
quantum simulators. 

In our opinion, the realization 

Let us consider the effective Hamiltonian obtained by the PAT scheme (20) for any dimensionality. We assume equal tun¬ 
nelings ta =: t, driving parameters Va =: r, and driving ratios r]a =• t], in all directions Va, and set (pa = 0. The limit 
of strong Hubbard repulsion in this case corresponds to tJriv) where the parity-violating tunnelings described by 

the terms and in Eqs. (63)-(64) can only take place virtually. As in the standard Hubbard model [100, 101], 
such virtual tunnelings can be calculated by a Schrieffer-Wolff-type unitary transformation ^eff = where S = 

—i(A's^d ~ ^d^s)/is responsible for eliminating the energetically forbidden tunnelings to second order in the small expan¬ 
sion parameter ^ = tJr (v) /^^U’ Considering the commutation properties of the different operators defined so far, one finds 
that Hsff = Hioc+Kq + Vm + ([A's^d,-Kd^s] + [^s^d,-Ko] + [^o,-*^d^s])/5t/|; + For hole-doping about half-filling, one 

can project onto the single-occupancy sub-band, which yields the aforementioned t-J model (66) with an additional density- 
dependent next-nearest-neighbor tunnelling 




^ = LL ^ (-^!a«i+e„.a/i+„“ <7 +^V4e„.a/i+e„.a/i+u“.cT + H-C.) 


oc,n 


(67) 


where we have introduced the effective cold-atom parameters 

f = rJo(r/), J=4t%{r]f/5Un, (68) 

and the next-nearest-neighbor vectors (e.g. for ID = 20^, for 2D {u^ = 2ex,U2 = = 203;, = 0^ + 63 ;}). Let 

us note that this additional tunnelling requires that the target site is populated with a hole, as otherwise Hence, 

close to half-filling (nt) ^ 1, this term is reduced by a factor (1 — (^/))/4 with respect to the Heisenberg super-exchange, and is 
typically neglected in the literature [92] . 
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As announced above, we have obtained an effective t-J model where the ratio of the coupling constants (68), namely J/t = 
4(^ 1^(7] )/Jo(t7 ), can be tuned by modifying the intensity of the moving optical lattice. Even if ^ < 0.1 in the regime of validity of 
the t-J model, we can make Jo(7]), such that exchange coupling is not required to be much smaller than the tunnelling 

as in the standard Fermi-Hubbard model. In this way, we can explore the full phase diagram of the t-J model. 

In the ID case, theoretical predictions about the phase diagram are supported by very accurate numerical methods [102]. 
Such results could serve to benchmark the accuracy of the proposed quantum simulator, which can be prepared in a regime 
corresponding to a metallic phase (i.e. a repulsive Luttinger liquid [103]), a gapless superconductor (i.e. an attractive Luttinger 
liquid [103]), or the so-called spin-gap phase [104], which consists of a quantum fluid of bound singlets with gapless density 
excitations, a gapped spin sector, and enhanced superconducting correlations (i.e. a Luther-Emery liquid [105]). Finally, for 
sufficiently strong super-exchange interactions, the antiferromagnetic order expels the doped holes leading to separated hole- 
rich and hole-poor regions (i.e. phase separation [106]). The richness of this phase diagram highlights the potential of our PAT 
scheme, and contrasts with the standard ID Fermi-Hubbard model where only the repulsive Futtinger liquid can be achieved. 
For instance, the Futher-Emery liquid, which has eluded experimental confirmation so far, requires 2.5? and thus lies out of 
the range of parameters that can be obtained from the repulsive Hubbard model. Moreover, our quantum simulator would allow 
to test the numerical results [107] predicting the disappearance of the phase separation in favor of enlarged superconducting and 
spin-gap regions, once the next-nearest-neighbor tunnelling terms in Eq. (67) are considered. 

In the 2D case, a detailed understanding of certain regions of the phase diagram is still an open problem, and the subject of 
considerable debate. As emphasized in [108], theoretical predictions are difficult to verify due to (/) the absence of controlled 
analytical methods, and (ii) the limitation of numerical methods to small system sizes where finite-size effects can affect the 
predicting power. From this perspective, the proposed quantum simulator may eventually address some of the following open 
questions regarding the properties of the t-J model. Variational methods based on a resonating-valence-bond trial state [16] (i.e. 
a linear superposition of all possible configurations of singlet pairs with a weight that depends on the pairing symmetry), have 
predicted a rich phase diagram [92] with regions off/) ferromagnetism, (ii) 5'-wave pairing, (Hi) J-wave pairing, (iv) coexistent 
antiferromagnetism and superconductivity, and fv) phase separation. However, this variational approach introduces a certain bias 
through the choice of the particular set of ansatzs, and this compromises its reliability leading to considerable controversy in the 
community [2]. In particular, there are contradictory predictions for low dopings and not too large ratios of///, which turns out 
to be the regime of interest for the high-Tc cuprates. For instance, the results of [92, 1 10] contradict those of [106, 109], regarding 
the onset of the phase separation in this low-doping region. There has also been some disagreement regarding the regions of 
superconductivity predicted by the variational approach [92], exact diagonalization [114], and quantum Monte Carlo [115]. From 
the perspective of the high-Tc cuprates, addressing the conflicting predictions in [111, 112] and [92, 113] about the existence 
of stripe phases (i.e. inhomogeneous charge and spin distributions) in the t-J model is even more compelling, as these have 
been measured experimentally in the cuprates. If the t-J model is to function as a canonical model of the cuprates, as advocated 
in [16, 92, 93], it is important to settle this dispute and determine if it admits stripe phases. We believe that the proposed 
cold-atom experiment could be helpful in this respect. 

Before closing this section, let us comment on two additional possibilities for the t-J model quantum simulator. The first, and 
most obvious one, is the possibility of controlling the spatial anisotropy of the parameters of the effective Hamiltonian (66) by 
simply exploiting the dependence of the dressed tunnelings on the lattice axes. Accordingly, the effective t-J model becomes 
anisotropic, such that the anisotropy of the tunnelings ta = /aJo(t]a), and the super-exchange couplings fa = 4/^J^(r]ct)^/5f/|^, 
can be controlled through the intensities of the static and moving optical lattices along the different axes. Such anisotropy 
becomes especially interesting in the context of certain cuprate ladder compounds [116], which can be modeled by a number 
i G {1,... of one-dimensional t-J chains that are coupled to each other by the transverse tunneling t' and super-exchange 
coupling /. This defines the so-called rungs of the t-J ladder Hamiltonian , where 


= 11 (-f'Acr/i+i.to +H.C.) +H/(5,-,r 

i,oi=l ^ ^ i i 

i,ai=l ^ ^ i i 


(69) 


Such ladder Hamiltonians can be implemented in our quantum simulator if we supplement the above scheme (67) with additional 
static lattices along the y-axis with commensurate wavelengths with respect to the original lattice. For instance, combining two 
lattices with doubled wavelengths Ay = 2Ay, leads to an array of decoupled ladders with = 2 legs [117]. By adding further 
harmonics, one may create ladders with other numbers of legs, at least in principle (e.g. the first n£ harmonics of the Fourier 
series of the square-wave function yield an approximation to an array of decoupled Tz^-legged ladders). Hence, the interaction- 
dependent PAT scheme leads to Eq. (69) with tunable number of legs and Hamiltonian parameters 


4t^]r{Tlx)^ 


4tpr{riy)^ 




t - tyJ()(Tly,), J 


(70) 
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together with the corresponding density-dependent next-to-nearest-neighbor tunneling (67), typically neglected for small dop¬ 
ings. These t-J ladders provide a very interesting interpolation between the well-understood ID t-J model, and the more 
intriguing 2D case full of open questions of relevance to high-Tc superconductivity. In fact, already for = 2 legs, the doped 
holes tend to pair [118] developing a superconducting d-wave-\\kQ order [119]. In addition to the phases that also occur for the 
ID t-J model [102], this J-wave superconductivity takes place in a wide region of the phase diagram [120], which includes the 
parameters relevant for the cuprates. Also, a very interesting even-odd effect reminiscent of the spin-gap presence/absence in the 
undoped system has been identified [121], whereby the hole J-wave pairing disappears for ladders with an odd numbers of legs. 
We note that some of these results [118] depend on ratios t'/J' < I that cannot be reached from standard ladder Hubbard models 
where d :^J'. Likewise, the regime / ^ and its connection to short-range resonating valence bond states [122, 123], can¬ 
not be reached from standard Hubbard models. It would be very interesting to test these predictions with our quantum simulator, 
which allows exploring all these parameter regimes. 

Let us move to the last possibility of the t-J model quantum simulator: inducing a Heisenberg-Ising anisotropy in the super¬ 
exchange interactions. This leads to a very interesting t-XXZ model described by the Hamiltonian Hfxxz = 


H,xxz = II +I +‘5^51;,J +/,(5?5? 

i,a <7 ^ i,a 


<z 

i+Ca 




(71) 


where the Ising J^ and fiip-fiop interaction strengths are generally different, such that the t-J model is recovered when J^=J^. 
To achieve such an effective model with our quantum simulator, we must employ the PAT by a state-dependent moving optical 
lattice, which yields the effective Hamiltonian (26) in ID, and (30) in 2D. We assume equal tunnelings ta =: L spin-dependent 
driving parameters % =: %, and ratios r]a^a =• ^o^ in all directions Va, and set (pcy,a = 0- One may observe in Fig. 5(f) that 

second-order super-exchange interactions will depend on the spin configuration. In fact, we find 


ta — Ho(t]<7), Jx 




(72) 


such that the Heisenberg-Ising anisotropy ^ = J_l/Jz can be tuned all the way from small spin quantum fiuctuations C ^ ^ 
(i.e. t-J^ model), to the isotropic regime where spin quantum fiuctuations play an important role ^ 1 (i.e. t-J model). These 

Hamiltonians have been studied in the context of the propagation of a single hole in an antiferromagnetic matrix. For instance, 
in the t-J^ model, the tunneling of the hole leaves behind a string of flipped Ising spins that costs an energy proportional to 
the string length, such that the holes are almost [124] localized to the site where they were doped [125]. The situation is 
considerably more complex as the spin fiuctuations are switched on ^ > 0, and the t-J model is approached C ^ 1 [126], 
and some controversy regarding the limitations of the different analytical or numerical methods has been identified [127]. The 
possibility of controlling the amount of spin fiuctuations in our quantum simulator, together with the possibility of using the 
high-resolution optics of quantum gas microscopes [128, 129] to create localized holes and watch them propagate in real time, 
opens a very nice perspective in accessing this quantum many-body effect with ultracold atoms in optical lattices. 


C. Synthetic dynamical Gauge fields 

Gauge theories play a prominent role in several areas of modern theoretical physics, such as the strong interactions between 
quarks and gluons in quantum chromodynamics. Although perturbative predictions are reliable at short distances, the low-energy 
regime eludes a perturbative treatment and leads to numerous unsolved questions such as the phase diagram of quark matter [3]. 
Therefore, a quantum simulator for quantum field theories of coupled Gauge and matter fields would indeed be very useful. 

From a general perspective, this would require (i) Gauge (matter) degrees of freedom evolving under the Hamiltonian field 
theory (H^), and (ii) a tunable interaction introduced by the so-called minimal coupling. So far, most of the theoret¬ 
ical and experimental progress has considered static/background Gauge fields = 0, where the Gauge is fixed [31]. There is 
however an increasing interest in promoting this situation to a regime of dynamical Gauge fields ^ 0, as reviewed in [130]. 
To preserve the Gauge symmetry, one parallels the construction of lattice Gauge theories [131]. In the Hamiltonian formula¬ 
tion [132], the fermionic matter field vr(r) ^ is defined on the sites of a J-dimensional lattice r ^ n = ia, where i is a 

vector of integers and a the lattice spacing, whereas the Gauge degrees of freedom are defined in terms of unitary matrices t/ij. 
Such unitaries can be expressed in terms of a Gauge field f/i j = defined on the links ry = (ri rj)/2 of two neighboring 

sites connected by ^ = (rj —ri)/a. The complete Hamiltonian for the lattice Gauge theory 

// = 7%-I?ij +H.C.)+IeiV/iVi, (73) 

(iJ> i 

requires a particular Gauge-invariant construction of and a particular choice of tunnelings and on-site energies £{, which 
lead to the corresponding Gauge and matter field theories in the continuum limit a 0 (e.g. minimally-coupled Maxwell and 
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Dirac quantum field theories). Moreover, one must enforce Gauss law by considering only the physical states of a sector of the 
Hilbert space. Although there are very interesting proposals to accomplish this goal in the pure Gauge sector [133], and the 
complete lattice Gauge theory [134], we will restrict to a simpler quantum simulator of Eq. (73) with a Hamiltonian that is 
not Gauge invariant. Although departing form the standard formulation of lattice Gauge theories, the dynamical character of the 
fields brought by ^ 0 can lead to very interesting phenomena, which might be still associated with a dynamical Gauge field 
theory within a fixed gauge [135]. Alternatively, one can use focus on lattice field theories that lead to interacting relativistic 
quantum field theories at low energies, such as the Thirring and Gross-Neveu models [136]. 


7. Interacting relativistic quantum field theories: Yukawa-type couplings 


Let us focus on the ID case, and consider the PAT of a Fermi gas by a spin-dependent moving lattice (26) in a parameter regime 
fulfilling Eq. (24) for % = 2, and SUyy = 0. Let us note that the same can be obtained for a Bose-Fermi mixture (54), provided 
that the hardcore constraint is considered. In order to obtain a quantum simulator of quantum matter coupled to dynamical 
gauge fields (73), one needs to find a particular set of parameters such that: (/) the tunnelling amplitude (27) becomes a simple 
c-number, and (//) the tunnelling phase in Eq. (26) is non-vanishing only for one of the pseudospin states. In order to fulfill (/), 
the ratio of the moving-lattice intensities with respect to the detuning must fulfill Jq (t]^) = J 2 which can be achieved with 
any of the following values rj^ G {4.89,8.29,11.53,...}. By direct substitution in Eq. (27), one finds that the dressed tunnelling 
strength becomes the desired c-number Ix^raAnj+i ^ = Ixh =• order to fulfill (//), it suffices to set the relative 

phase of the moving optical lattice (p; = 0 for the pseudospins that will play the role of the Gauge field. Hence, the tunnelling 
phases in Eq. (26) fulfill 9 ^ = 0, and 2(py =: —0 7 ^ 0, such that the atoms with pseudospin ||) play the role of the dynamical 
Gauge field for the atoms with pseudospin |t) (see Fig. 7), namely 
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This is the general result of this section. The PAT scheme has allowed us to build an effective Hamiltonian where the atoms with 
one of the pseudospins hop freely in the lattice and play the role of a dynamical ’Gauge’ field for the atoms with the remaining 
pseudospin. This must be contrasted with other interesting proposals [51, 54], where the tunnelling Peierls phase for atoms 
depends on their own density, such that the roles of matter and Gauge fields cannot be distinguished. We believe that these type 
of Hamiltonians (74), and their straightforward generalization to higher dimensions, to other lattices, to hardcore bosons, or to 
slightly modified PAT schemes [137], will lead to several interesting many-body phenomena that deserve to be studied in further 
detail. To illustrate this richness, we describe a particular example that leads to an interesting relativistic quantum field theory in 
the continuum. 

The effective Hamiltonian (74) corresponds exactly to the structure of the dynamical gauge field theory in Eq. (73) if we 
make the following identifications: (i) The fermionic quantum matter is represented by the atoms with pseudospin |t), namely 
¥i •= Gauge degrees of freedom will be some collective low-energy excitations of the atoms with pseudospin ||). 

In the half-filled ID case, it is well-known that the particle-hole excitations of the free tight-binding Hamiltonian can be mapped 
onto a pair of bosonic branches [103, 138], which will play the role of the gauge degrees of freedom 




Y^b^ ^ -|~ bq Y/bq L / 5 


(75) 


where we have introduce the quasi momentum q = Inn/L for n £ Z+, the effective speed of light and the 

bosonic operators for the particle-hole excitations {bqy) around the right (left) Fermi point. Moreover, we have assumed 
that the modifications with respect to half-filling coming from the weak parabolic trapping, encoded in | 
can be accounted for using a local chemical potential 5/7/. Note that, in the continuum limit, this bosonic Hamiltonian becomes 
a scalar field theory with the energy zero set at namely 

where the scalar field 0(v) = 0 r(x) — 0l(x) is expressed in terms of the inverse Fourier transform of the bosonic operators bq^R, 
bq^Y and their Hermitian adjoints, and 7l{x) is its canonically-conjugate momentum. Although there is no Gauge invariance in 
the massless scalar field theory (76), the dynamical bosonic field will interact with the fermionic matter, and can still lead to 
interesting phenomena. The particular form of the interaction brings us to the final identification (in) the gauge field that dresses 
the tunnelling of the quantum matter in Eq. (73) becomes aAx{xiq^i) = 0(^/+i ^ which amounts to the density difference 

of the ’Gauge’ species. In order to express the gauge unitaries in terms of the bosonic particle-hole excitations, we use 


aAx{xiq^i) =ae 


^^•+14 


= a^Gdx {x)\i/{x)) -^Ajc{x) = 6ad^(j)(x), Ui^^i = 1 -fiiaAjc{x) + ^(a^A^), 


(77) 


a 



26 


where we have again applied the continuum limit. To be consistent with such a limit, we note that the bare tunnelling of the 
fermionic quantum matter corresponds to the ID version [139] of the so-called Kogut-Susskind fermions [132]. At half-filling, 
one obtains the massless Dirac field theory in the continuum limit, which minimally couples to a derivative of the scalar field 

H^qfy = H^ + f ^(^'I'[(x)(5 + ct(-ia,-A,(x)))4'L(x)+'P^(x)(5 + ct(+i5,+A,(x)))'PR(x)^, (78) 

where we have introduced the fermionic field operators T^r(v) (T^l(-^)) for the right (left) moving fermions with pseudospin 
(7 =t, the energy difference due to the different Zeeman shifts 5 = 8 ^ — and the effective speed of light We have 

obtained a peculiar relativistic theory of interacting quantum fields. Instead of the standard Yukawa coupling between scalar 
and fermionic fields 1/701//, we get a minimal coupling with the derivative of the scalar field 1/77^5^01//. Moreover, the effective 
speeds of light of the scalar and fermionic particles can be tuned independently. Anyhow, scattering between the fermions 
will occur due to the exchange of scalar particles, such that the cold-atom experiment could be exploited to calculate scattering 
amplitudes in the spirit of [140]. However, it is not clear how one would create the initial incoming particles and measure the 
outgoing scattering probabilities in our scheme. A simpler goal would be to study collective properties of the model (78). For 
instance, for very large 5, fermion-fermion interactions will be mediated by the virtual exchange of scalar particles, such that 
the fermionic properties of the groundstate will be modified (e.g. correlation functions). Increasing the fiux 6 could lead to new 
phases departing from the Luttinger-liquid phase of the free Kogut-Susskind fermions. All these questions could be studied with 
the proposed quantum simulator. 


2. Correlated topological insulators: Hofstadter-type models with dynamical Gauge fields 


Topological insulators represent a family of holographic phases of matter that are insulating in the bulk and conducting at 
the boundaries [141]. These states are topologically different from trivial band insulators, as the bulk is characterized by a 
finite topological invariant that cannot be changed unless the bulk energy gap is closed (i.e. phase transition). Another difference 
occurs at the boundaries, where gapless edge states are responsible for the conductance. In the absence of a boundary energy gap, 
it is the chirality, or some additional symmetry of the problem, which underlies the robustness of the conductance. This is clearly 
exemplified by the so-called Hofstadter model [142], which describes fermions in a square lattice subjected to a perpendicular 
magnetic field, and displays the aforementioned bulk [143] and edge [144] properties. Similar phenomenology also arises in 
the Haldane model [145], which is a topological insulator in the same symmetry class (i.e. time-reversal symmetry breaking of 
the integer quantum Hall effect). Remarkably, other instances of topological insulators belonging to different symmetry classes 
have also been found, such as the Kane-Mele model [146] built from two time-reversed copies of the Haldane model, or the 
time-reversal Hofstadter model [147] built from two time-reversed copies of the Hofstadter model. 

Paralleling the effect of interactions in the quantum Hall effect [148], one expects that even more exotic phases of matter 
will appear when considering the effect of correlations in the above models [149]. So far, the typical route to introduce such 
correlation effects has been to include the effect of on-site and nearest-neighbor Hubbard interactions in the Haldane [150], 
Kane-Mele [151], or time-reversal Hofstadter [152] models. These studies show that the topological features are robust to 
interactions, but no other exotic phases such as topological Mott insulators [153] (i.e. interaction-induced bulk gap and protected 
edge states) or topological fractional insulators [154] (i.e. fractional excitations and protected edge states) were found. In this 
section, we consider explicitly the 2D PAT by a spin-dependent moving lattice (28), and discuss how this scheme may be used to 
explore a new type of correlation effects introduced by substituting the fixed background gauge field by a dynamical one in the 
standard [142], and time-reversal [147], Hofstadter models. Given the recent realizations of both models with non-interacting 
cold atoms [44], we believe that future experiments will be able to explore the full phase diagram, and the possibility of finding 
more exotic phases brought by the interactions with the dynamical Gauge field. 

(i) Hofstadter model in a dynamical Gauge field - Let us consider the 2D scheme (28) leading to the effective Hamilto¬ 
nian (30). We generalize our prescription for the ID case, and set for = 2, and 5Uyi = 0. In order to make the dressed 
tunnelling amplitude a c-number, we tune again the ratio of the moving-lattice intensities with respect to the detunings to 
Vx — Vy = ^ {4.89,8.29,11.53,...}, such that Jq (t]^) = J2 We can thus define the effective tunnelling amplitudes 

along the v- and y-axes as t^^^^ := 4 J 0 and := tyJo Additionally, we need to control the relative phases of 

the moving lattices, such that only one of the pseudospins develops a non-vanishing Peierls phase = 0, but ^ 0. 

Accordingly, the effective Hamiltonian (30) becomes a Hofstadter model in the Landau gauge for the atoms with pseudospin |t) 


F/eff = ^ 






(79) 


Here, the synthetic Gauge field according to Eq. (73) corresponds to aA{x[^i^cJ = 0(^i+e;„; — and thus depends on the 
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atoms with the remaining pseudospin ||), which evolve under the free tight-binding Hamiltonian 


= -I (4^4/i+e.4+4.4/i+e,4+H-C-) • 


(80) 


The same can be obtained for a Bose-Fermi mixture (54), provided that the hardcore constraint is considered, such that the 
hardcore bosons play the role of the Gauge field. Let us note that the dynamical Peierls phase cannot be gauged away due to 
its inhomogeneity, and thus corresponds to a non-trivial ’Gauge’ field. Alternatively, one can compute the Wilson loop operator 
around a square plaquette. 


W(^ — f/i+ey,i^i+ex+e3;,i+e3;^i+e;c,i+ex+ey^i,i+ex 


— ^i+ex,i ^i+ey,t+^i+ex+e3;,|) _ 
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and check that it is a non-trivial operator for 6 G (0,2;r). As the ’Gauge’ fields commute at any position, one can regard the 
Hamiltonian (79) as a dynamical Abelian Hofstadter model. The quantum simulator will be able to explore this interesting 
model, and the fate of the Hofstadter quantum Hall phase, in the light of the aforementioned correlated topological insulators. 

(ii) Time-reversal Hofstadter model in a dynamical gauge field - We can now generalize the above construction to a time- 
reversal invariant situation, such as having two copies of the Hofstadter model subjected to anti-parallel magnetic fields. In 
our case, this is straightforward if we tune the relative phases as follows = 20 7 ^ 0. Moreover, we also set 

^eff a ~ ^eff a ^eff,a by controlling the moving optical lattices appropriately. Accordingly, the effective Hamiltonian (30) 
becomes a doubled Hofstadter model in the Landau gauge 
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where atoms of any pseudospin play the role of the dynamical ’Gauge’ field for the atoms with the remaining pseudospin, and 
the magnetic fluxes are opposite for each pseudospin. Accordingly, the time-reversal symmetry should invert the flux 6^—6, 
and flip the pseudospin /i Hence, the quantum simulator can explore the fate of the fate of the time-reversal Hofstadter 

quantum spin-Hall phase due to the presence of strong correlations, and dynamical effects of the Gauge field. 


IV. CONCLUSIONS AND OUTLOOK 

We have introduced an interaction-dependent photon-assisted tunneling by combining the strong Hubbard interactions of cold 
atoms in optical lattices with a periodic driving stemming from a moving optical lattice. This effect leads to exotic Bose-, Fermi-, 
and Bose-Fermi Hubbard models with generalized tunnelings whose strength depends on the atomic density through a Bessel 
function that is controlled by the intensity of the moving lattice. Additionally, the effective Peierls phase of the tunnelling also 
depends on the atomic density, but through the phase of the moving lattice. We have argued that this effect can be exploited as a 
flexible tool to implement a variety of quantum simulations of quantum many-body models in the context of strongly-correlated 
electrons and high-energy physics. In particular, our scheme may allow to explore paradigmatic models, such as the t-J model, 
in regimes that were previously inaccessible to cold-atom experiments. Such an experiment would be very relevant to test the 
accuracy of current approximate methods that study the phase diagram of the model, and the possibility of displaying d-wave 
superconductivity mediated by strong correlations. Moreover, this proposal introduces a new perspective in the realization of 
dynamical Gauge fields departing from the lattice Gauge theory approach, which can also lead to interesting quantum many- 
body models. We once again remark that such the dynamics of such synthetic Gauge fields is not itself Gauge invariant, and thus 
such models cannot be related to a lattice Gauge theory. 

From an experimental point of view, the quantum simulation of the models that does not require super-exchange (e.g. phases 
related to the bond-charge interactions, and Nagaoka ferromagnetism), will require less stringent time-scales as the dressed 
tunneling can be made of the same order of magnitude as the bare one by choosing the right parameters. Conversely, the quantum 
simulation of the t-J model will require much slower timescales of the adiabatic protocol, and also stringent cooling conditions 
to guarantee that the adiabatic protocol is closer to idealized one. However, the dressed super-exchange (70) can be made also 
on the same order of magnitude of the bare super-exchange if the right parameters are chosen. For the dynamical Gauge fields, 
although the timescales would be favorable as one is interested in the dressed tunneling and not in second-order processes, the 
scheme gets complicated by the requirements of the photon-assisted-tunneling scheme (e.g. state-dependent moving lattices). 

Although we have focused on cold atoms, the scheme can be applied to other setups, provided that the relevant dynamics can 
be described by a lattice model, and one can control the periodic driving and the strong interactions. This might be the case of 
trapped-ion crystals, where the local vibrations and electronic states lead to a lattice model with bosonic and pseudospin degrees 
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of freedom, and the interactions and periodic drivings are provided by their interaction with laser beams.A similar situation 
arises for superconducting circuits by considering photons in arrays of microwave cavities and superconducting qubits leading 
to the aforementioned lattice models. 
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